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1.  Executive  Summary 

1.1  Introduction 


This  report  describes  work  performed  from  1  April  2004  to  30  May  2008  under  AFOSR  grant 
number  FA9550-04-1-0219  entitled  “Optical  Pulse  Interactions  in  Nonlinear  Excited  State 
Materials”.  We  cover  five  general  areas:  (1)  optical  limiting,  (2)  organic  multiphoton  absorber 
(AFX),  (3)  numerical  methods,  (4)  relationship  between  the  building  blocks/decay  rate  matrix 
and  the  coupled  propagation/rate  equations,  and  (5)  semiconductor  quantum  dots.  The 
introduction  to  the  project  and  an  overview  of  each  area  is  given  in  Section  2.  Then  in  Section  3 
we  present  the  results  for  each  area.  Section  4  contains  the  conclusion  and  Section  5  lists  the 
references.  There  are  three  detailed  appendices  of  the  mathematics  at  the  end. 

1.2  List  of  People  involved 

M.  Potasek,  PI 

D.  McLaughlin,  Co-PI 

E.  Parilov,  research  associate 
Other: 

Consultative  and  advisory  functions  to  other  laboratories  and  agencies 

•  Consultation  and  advisory  functions  to  Dr.  Chris  Brewer,  WPAFB/ML, 

Dr.  Ruth  Pachter,  WPAFB/ML,  Dr.  Paul  Fleitz,  WPAFB/ML,  Dr.  Mark 
Walker 

•  Transfer  of  documentation  and  results  on  multiple  pulses  in  femtosecond 
microlithography,  characterization  of  photosensitive  materials  for 
inducing  and  fabricating  novel  devices  such  as  photonic  crystals,  sensors, 
biosensors,  waveguides,  MEMS,  including  numerous  conference  calls  and 
e-mail  communications  discussing  results  with  researchers  at  WPAFB 

Transitions 

•  Transition  of  new  findings  and  all  related  materials  including  related 
publications  to  WPAF/ML. 

1.3  List  of  Publications/Reports/Presentations 

(a)  Referred  Publications: 

E.  Parilov  and  M.  Potasek, "Generalized  theoretical  treatment  and  numerical  method  of  time- 
resolved  radially  dependent  laser  pulses  interacting  with  multiphoton  absorbers,"  J.  Opt.  Soc 
Amer.  B  23,  1894  (2006), 

E.  Parilov  and  M.  Potasek, "Generalized  theoretical  treatment  and  numerical  method  of  time- 
resolved  radially  dependent  laser  pulses  interacting  with  multiphoton  absorbers,"  Virtual  Journal 
of  Ultrafast  Science,  Oct.  2006. 

(b)  Reports: 


1 


M.  Potasek  and  D.  McLaughlin,  “Optical  Pulse  Interactions  in  Nonlinear  Excited  State 
Materials”,  AFOSR  Reports,  2004-2008. 

E.  Parilov  and  M.  Potasek,  “Unifying  Approach  to  Numerical  Method  for  an  Electromagnetic 
Wave  Propagating  Through  a  Photoactive  Material”,  Aug.  2007. 

M.  Potasek  and  E.  Parilov,  “Semiconductor  Quantum  Dots:  Experiments  and  Theory”,  May  30, 
2008. 

(c)  Presentations: 

M.  Potasek  and  E.  Parilov,  “Generalized  theoretical  treatment  and  numerical  method  of  time- 
resolved  radially-dependent  laser  pulses  interacting  with  multi-photon  absorbers”,  Univ.  of 
Arizona,  Oct.  17,  2006. 

(d)  Participation  in  Conferences: 

•  International  Symposium  on  Optical  Science  and  Technology,  49th  Annual 
Meeting,  2-6  August  2004,  Denver,  Colorado 

•  Conference  on  Two  Photon  Absorbers  and  Nonlinear  Optics,  Washington, 
DC,  Jan  20,  2004 

•  Photonics  West,  San  Jose  Convention  Center,  San  Jose,  CA,  Jan.  21-26, 
2006 

•  Photonics  West,  BIOS,  LASE,  San  Jose  Convention  Center,  San  Jose,  CA, 
Jan.  23-25,  2007 

•  Photonics  West,  BIOS,  OPTO,  San  Jose  Convention  Center,  San  Jose, 
CA,  Jan.  22-26,  2008 

1.4  List  of  Patents 

E.  Parilov  and  M.  Potasek,  “Method,  System  and  Software  Arrangement  for  Detennining  an 
Interaction  Between  and  an  Electromagnetic  Radiation  and  a  Material”,  US2007/0290147  Al, 
Dec.  2007. 

2.  Introduction  to  the  Project 
2.1  Overview  of  Optical  Limiting 

The  general  area  of  nonlinear  absorbers/optical  limiting  has  grown  substantially  over  the 
last  two  decades.  There  is  increasing  interest  in  various  mechanisms  of  optical  limiting.  A 
number  of  review  articles  and  conference  reports  demonstrate  the  vitality  of  this  subject1'11.  In 
general,  an  optical  limiter  keeps  the  power,  intensity,  energy  or  energy  density  transmitted  by  an 
optical  system  below  a  predetermined  maximum  value  that  is  independent  of  the  size  of  the  input 
pulse  while  maintaining  a  high  transmittance  at  low  input  power.  The  many  applications  of  the 
device  include  laser  power  regulation,  laser  mode-locking,  optical  pulse  shaping,  signal  level 
processing,  and  sensor/detector  protection. 

Importantly,  the  device  must  operate  within  a  fast  enough  response  time  to  protect  a  sensor 
or  regulate  pulsed  laser  sources.  The  transmission  must  be  reduced  in  a  time  period  shorter  than 
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the  pulse  width.  With  laser  pulse  durations  ranging  from  milliseconds  to  femtoseconds, 
electrooptic  devices  such  as  electrooptic  switches  that  operate  on  the  nanosecond  time  scale  are 
often  too  slow.  On  the  other  hand,  passive  devices  use  a  nonlinear  optical  medium  that  has  a  fast 
response  time  and  acts  as  a  sensor,  processor  and  modulator  .  However,  typical  nonlinear 
materials  are  not  sufficient  optical  limiters  by  themselves.  As  a  result,  the  passive  system  places 
the  nonlinear  optical  components  in  or  near  a  focal  plane.  The  focal  point  should  be  located  near 
the  exit  plane  of  the  material  to  utilize  the  entire  path  length.  The  time  scale  for  various  optical 
limiting  methods  is  shown  in  Fig.  1. 
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Fig.  1  Time  scales  for  various  optical  limiting  effects. 

Optical  limiting  now  requires  a  sophisticated  set  of  properties  including  various  nonlinear 
materials,  different  physical  processes,  various  device  designs  including  a  number  of  lenses  and 
other  optical  devices,  optical  pulse  shaping,  and  the  characteristics  of  the  laser  source  including 
wavelength,  pulse  duration  and  pulse  repetition  rate.  This  range  of  materials,  processes,  and 
devices  involve  complex  optical  propagation  requiring  a  vast  diversity  of  parameters  that  cannot 
be  fully  analyzed  without  the  use  of  various  numerical  codes.  Along  with  new  materials  and 
device  structures,  new  numerical  codes  will  be  key  elements  in  the  development  of  optical 
limiters.  Numerical  investigation  will  be  required  to  make  additional  advances  by  examining  the 
various  combinations  of  material,  physical,  and  device  parameters  that  would  be  too  time 
consuming  and/or  costly  to  perform  in  the  laboratory. 

Physical  mechanisms  used  in  optical  limiters  include  nonlinear  refraction,  nonlinear 
scattering,  two  photon  absorption,  thennal  nonlinearity,  and  reverse  saturable  absorption. 
Materials  used  in  optical  limiters  include  fullerenes,  porphyrins,  liquid  crystals,  carbon  black  and 
nanostructured  materials1'11.  Various  device  design  and  structures  also  play  an  important  role. 
Examples  of  different  device  design  include  tandem  and  graded  density  structures7. 

2.2  Overview  of  Organic  Multiphoton  Absorbers  (AFX) 

The  number  of  applications  for  materials  exhibiting  nonlinear  multiphoton  behavior  has 

12 

grown  significantly  in  the  last  two  decades.  They  include  safe,  high  precision  diagnostic  tools, 
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effective  treatments  for  various  cancer  diseases,13  optical  limiters,14'18  biological  detectors,19 
three-dimensional  microfabrication,"  '  fluorescent  imaging  systems,"  and  optical  storage."  ’ 
Because  of  their  importance  for  a  wide  range  of  applications,  the  design  of  molecules  with 
large  two-photon  absorption  (TP A)  is  highly  significant.26  Electron  donor-acceptor  systems  with 

Tf-conjugation  are  often  used  to  enhance  the  TPA. 26,27  However,  it  was  recognized  in  the  1970s 

that  the  TPA  is  often  followed  by  excited  state  absorption  (ESA)."  Furthermore,  nonlinear 
transmission  measurements  and  Z-scan  measurements  of  organic  materials  show  the  presence  of 

ESA.  ’  Recently,  a  family  of  D-ir-A  molecules  has  been  introduced  for  TPA  applications, 

known  as  the  AFX  .  These  chromophores  are  gathering  significant  interest  for  applications  in 
multi-photon  processes. 

For  the  AFX  group  of  nonlinear  chromophores,  ’  "  such  simulation  is  based  on  the  solution 
of  the  coupled  system  of  propagation  and  rate  equations  with  the  latter  described  by  a 
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phenomenological  five-level  absorption  model  depicted  in  Fig.2. 
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Fig.  2.  Energy-level  diagram  for  AFX  chromophores.  Electronic  states  are  labeled 
N0  -  N4 .  The  upward  arrows  correspond  to  photoexcitation  transitions  and  downward  arrows 

correspond  to  electron  decays.  The  cr’s  are  absorption  cross-sections  and  the  k  ’s  are  the  decay 
rates. 


Propagated  light  may  decrease  due  to  TPA,  which  excites  electrons  from  the  ground  level 
N0  to  the  singlet  excited  state  A, ,  and/or  due  to  ESA  from  singlet  and/or  triplets  states  Nl,N1. 

Depending  on  pulse  width  and  intensity  of  the  incident  light,  the  electron  population  densities 
change  which,  in  turn,  changes  the  transmittance  characteristics  of  the  material. 

2.3  Overview  of  Numerical  Method 

To  analyze  and  predict  the  experimentally  observable  behavior  of  such  materials  under  laser 
irradiation,  one  may  need  to  solve  a  nonlinear  system  of  differential  equations.  Because  not 
every  system  may  be  reduced  to  a  simple  fonn  and  be  solved  analytically,  it  is  extremely 
important  to  have  effective  and  robust  numerical  simulation  tools.  Approximate  theories  of 
propagation/transmission  through  nonlinear  materials  were  designed  by  researchers  in 
conjunction  with  their  particular  experimental  data.  Many  of  these  approximate  theories 
required  numerical  solutions;  therefore,  in  order  to  reduce  the  computation  time  several 
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researchers  have  used  approximate  analytic  expressions  based  on  some  simplifying  assumptions. 
Additionally,  because  traditional  approaches  are  often  used  in  conjugation  with  specific  laser 
systems  (e.g.  specific  wavelength  and  pulse  duration)  the  theoretical/numerical  analysis  has 
limited  applicability.  This  aspect  often  limits  the  predictive  capabilities  of  the  analysis  (e.g.  a 
theoretical  description  for  a  ns  pulsed  laser  cannot  predict  the  outcome  for  a  ps  or  fs  duration 
laser  system).  While  traditional  approaches  have  lead  to  agreement  with  specific  experiments  for 
specific  materials,  and  have  yielded  significant  insight,  particularly  in  absorbers  described  by 
single  photons  processes,  we  showed  that  these  approaches  must  be  refined  and  expanded  to 
meet  the  evolving  field  of  laser  transmission  through  excited  state  absorbers.  In  particular,  we 
provided  a  more  general  approach  for  calculating  the  transmission  of  laser  light  through 
nonlinear  absorbers.  We  use  measured  parameters  -  the  absorption  cross  sections  and  decay 
rates  -  in  our  analysis  and  do  not  attempt  to  derive  these  parameters  either  from  microscopic 
quantum  mechanical  theories  or  from  experimentally  measured  transmittance  data.  Rather  we 
establish  and  derive  algorithms  for  the  numerical  calculation  of  macroscopic  rate  equations 
including  the  assumptions  and  range  of  applicability  of  our  approach.  In  order  to  achieve  a  more 
comprehensive  numerical  method,  we  introduce  the  concept  of  “basic  (computational)  building 
block  diagrams”. 

We  described  a  general  numerical  method,  done  in  the  spirit  of  the  numerical  study  by 
Potasek  et  al.,16  to  evaluate  the  propagation  of  a  short-pulsed  laser  beam  in  a  multi-level  multi¬ 
photon  absorbing  material  that  includes  both  multi-photon  absorption  and  one  or  more  single 
photon  excited  state  absorptions.  We  outlined  a  broad  spectrum  of  non-linear  materials  to  which 
our  method  can  be  readily  applied,  given  the  photophysical  properties  of  the  materials.  Our 
calculations  can  be  perfonned  for  pulse  widths,  which  extend  from  the  ns  to  the  fs  regime.  We 
validated  our  technique  on  a  series  of  nonlinear  materials.  Because  our  method  is  applicable  to  a 
wide  range  of  materials  and  laser  systems,  it  may  enable  one  to  use  calculations  in  place  of  some 
experiments,  which  are  costly  to  perform  in  the  laboratory.  A  new  feature  of  our  numerical  tool 
enables  us  to  determine  exact  contributions  of  each  electronic  level  to  the  total  absorption,  which 
makes  it  possible  to  better  explain  the  roles  and  relationships  between  them  in  complex  multi¬ 
photon  absorbers. 

2.4  Overview  of  the  Relationship  between  the  Building  Blocks/Decay  Rate  Matrix  and  the 
Coupled  Propagation/Rate  Equations 

A  more  general  formulation  of  the  computational  “ building  blocks ”  shown  in  Fig.  3  can  be  used 
to  provide  information  and/or  data  relating  to  interactions  between  electromagnetic  waves  and 
“generic”  photoactive  materials.  For  example,  an  “ energy  lever  can  refer  to  one  or  a  group  of 
energy  states  which  can  contribute  to  macroscopic  physical  phenomena  such  as  those  described 
here  using  mathematical  models.  Energy  levels  may  apply  to  atoms,  molecules,  and/or  solids. 
Such  energy  levels  can  include,  but  are  not  limited  to,  electronic,  vibrational,  rotational,  and 
continuum  levels  of  an  atom,  molecule,  and/or  solid,  a  conduction  band,  and  a  valence  band  of 
semiconductors  and  metals.  The  computational  “building  blocks”  refer  to  an  absorption  or 
inverse  absorption  (stimulated  emission)  of  an  entity  including,  but  not  limited  to,  an  electron  or 
an  exciton,  between  different  energy  levels.  The  corresponding  decay  rates  appear  in  the  matrix 
described  below. 

2.5  Overview  of  Semiconductor  Quantum  Dot 

A  quantum  dot  (QD)  is  a  quasi-zero-dimensional  object  where  the  carrier  movement  is  restricted 
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in  three  dimensions.  The  bulk  crystalline  structure  of  the  semiconductor  is  maintained  in  the 
QD,  but  they  have  atomic-like  discrete  energy  states  due  to  3D  confinement.  Carrier  dynamics  in 
QDs  are  quite  different  from  bulk  materials  due  to  the  large  energy  level  spacing  and  enhanced 
surface-to-volume  ratio. 

There  are  several  generation/recombination  processes  in  semiconductors  which  are  described 
below.  A  particular  process  known  as  the  Auger  process  is  non-radiative  and  is  mediated  by 
electron-electron  (e-e)  Coulomb  interactions,  which  differ  amongst  atoms,  bulk  semiconductors 
and  QDs.  In  atomic  systems  the  e-e  coupling  is  much  stronger  than  the  electron-phonon 
coupling  which  makes  the  rates  of  Auger  transitions  much  greater  than  the  radiative  transitions. 
Thus  the  decay  of  multi-electron  levels  is  dominated  by  Auger  processes.  However,  in  bulk 
semiconductors  the  efficiency  of  the  Auger  process  is  reduced  because  the  interactions  of  the 
electrons  are  smaller  than  for  QDs.  The  situation  in  QDs  differs  from  the  bulk  in  that  the 
quantum  confinement  enhances  the  e-e  leading  to  increased  Auger  rates.  Additionally,  QDs 
have  a  large  surface  to  volume  ratio  that  leads  to  trap  states  in  the  band  gap.  In  many  QDs  the 
surface  states  are  passivated  by  capping  the  QD  with  a  different  semiconductor  material  or  an 
organic  material. 

3.  Results 

3.1  General  Numerical  Method 

3.1.1  Derivation  of  the  theory 

The  wave  equation  in  the  presence  of  the  electric  polarization  vector  is  given  by  Maxwell’s 
equation  in  the  scalar  form 


V2£(z,  r,  t)-\^E(z, r,  t)  =  1 2  f  2  P(z,  r,  t) . 
c  dt  £0<r  dt 

where  we  assume  that  V  •  E  =  0  ,  s0  is  the  pennittivity  and  c  is  the  speed  of  light  in  vacuum, 
optical  field  E  and  the  induced  nonlinear  polarization  P  are  given  by 


(1) 

The 


Ec  (z,  r,  t )  =  ^  E(z,  r,  t )  exp  [-/  -  k0z)]  +  c.c. 

Pc  ( z,r,t)  =  P{z,  r,  t )  exp  [-/  ( co0t  -  k0z)]  +  c.c. 


(2) 


where  co0(k0)  is  the  frequency  (wave  number)  of  the  incident  electromagnetic  field,  and 
E(z,r,t )  and  P(z,r,t)  are  the  slowly  varying  envelope  of  the  electromagnetic  field  and  the 
slowly  varying  polarization  vector,  respectively.  Then  by  using  the  well-know  slowly  varying 
envelope  approximation  (SVEA),  i.e.  the  pulse  envelope  varies  slowly  in  time  compared  to  an 
optical  period,  and  the  paraxial  approximation,33  i.e.  assuming  that  the  envelope  is  slowly 
varying  along  the  propagation  direction,  we  can  rewrite  the  scalar  wave  equation  in  the  form 


|A+is_^v; 

V  dz  c  dt  2k 


E(z,r,t)  =  —P(z,r,t), 


(3) 


where  V,  is  the  operator  for  the  transverse  variables  and  k  is  the  frequency  dependent  wave 


vector.  The  intensity  of  the  light  is  defined  by  I(z,r,t)  =  cns0  E{z,r,t ) 


where  n  is  the  linear 
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index  of  refraction,  and  the  photon  flux  at  the  carrier  frequency  co0  is  defined  by 
0(z,r,t)  =  l(z,r,t)/ha)o.  The  incident  intensity  of  the  laser  pulse  is  given  by 
I(z  =  0,r,t)  =  I0f(r,t )  where  f(r,t)  describes  the  radial  and  temporal  shape  of  the  incident 
pulse  whose  peak  value  is  I0 . 

The  polarization  vector  for  the  chromophore  can  be  described  by  the  density  matrix.34"37 
The  Hamiltonian  H  of  the  system  is  described  by  the  unperturbed  Hamiltonian  H0  and  the 

Hamiltonian  Hmt  (i.e.  perturbation  term)  such  that  H  =  H0  +  Hmt .  Hint  expresses  the  interaction 
of  light  with  the  molecular  system  by  using  the  electric-dipole  approximation7 
Hmt  =  -dxEx  =-d  E ,  where  d  is  the  electric-dipole  operator,  assuming  that  a  laser  is  linearly 
polarized,  and  d  and  E  both  are  aligned  along  x  axis.  The  eigenvalue  and  eigenenergy  for  the 
n-th  level  is  |  n)  and  en ,  respectively,  then  the  wave  function  and  unperturbed  Hamiltonian  can 

be  written  as  \y/)  =  \n),  H a\n)  =  en\n) .  The  density  matrix  operator  is  defined  as 

n 

g  =  \  ys){y/ 1  ?  and  the  equation  of  motion  is 


dg  ~iTl y  dgji 

dt  h  dt 


jk 


fl 


Z(HH  Si.-g.H,,). 


(4) 


where  the  matrix  elements  gik  correspond  to  a  polarization  induced  by  a  transition  between 


energy  levels  j  and  k  .  Photons  from  the  incident  light  are  absorbed  by  the  molecular  system 
which  enables  the  atoms/molecules  to  occupy  their  excited  energy  states  (electron,  vibrational, 
and  rotational).  Then  the  polarization  for  na  atoms/molecules  per  unit  volume,  where  R 
describes  the  distance  of  the  separated  charges  in  the  dipole  moment  approximation,  is  given  by 


P  =  na  e •  J yy  RdR  =  naJ^gjkdkj  =  naTr(dg)  ,  (5) 

jk 

where  e  is  the  unit  electric  charge.  After  excitation  to  higher  energy  states,  the  system  relaxes  to 
the  ground  state  through  radiative  and  nonradiative  transitions.  In  the  semi-classical  case  the 
relaxation  terms  are  added  to  the  equations  of  motion  of  the  density  matrix  phenomenologically. 
In  general,  the  quantum  mechanical  calculations  of  the  relaxation  terms  are  quite  involved.39'41 

For  an  understanding  of  the  behavior  of  an  ensemble  of  radiators  (e.g.  atoms,  molecules, 
excitons,  impurities  in  a  crystal)  in  a  field  of  resonant  or  nonresonant  electromagnetic  waves  it  is 
beneficial  to  characterize  the  time  scales  of  the  various  processes.  The  interaction  of  the  radiators 
with  non-resonant  atoms  (solvent)  can  give  rise  to  a  dephasing  rate,  yjk : 


d§jk 

dt 


-(rll,+i«>Jl)gJl  +A(»?s.  ~g,K)- 

M  l 


(6) 


3,1.2  Description  of  the  Generalized  Model 

As  laser  intensities  increase  many  organic  molecules  exhibit  multiphoton  absorption  involving 
both  singlet  and  triplet  states,  which  are  very  challenging  for  an  approach  based  solely  on 
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quantum  calculations.  As  a  result  a  phenomenological  model  based  on  spectroscopic  and  kinetic 
data,  has  lead  to  a  description  of  nonlinear  absorbers  involving  state  diagrams  or  Jablonski 
diagrams.42  This  model  provides  a  representation  for  population  dynamics  leading  to  the  so- 
called  rate  equations.  The  experimental  data  used  in  this  approach  includes  the  cross  section  of 
absorption  and  decay  rates  of  the  various  electronic  levels. 

Using  the  analysis  described  so  far  allows  us  to  encompass  a  large  group  of  nonlinear 
materials  which  may  be  characterized  by  applying  one  general  numerical  method.  In  fact,  there 
are  several  types  of  absorption  mechanisms  within  nonlinear  materials  depending  on  the  number 
of  photons  absorbed  simultaneously  and  on  which  states  absorption  occurs.  Even  though  our 
method  can  be  used  for  N-photon  absorbers  with  both  singlet  and  triplet  levels,  we  limit  our 
examples  to  at  most  three  photon  absorption  materials  in  order  to  compare  our  numerical  method 
to  experiments.  In  this  way  we  can  distinguish  five  different  types  of  absorption  mechanisms 
which  are  depicted  on  Fig.  3(B0)-(B4)  as  transition  diagrams.  According  to  the  selection  rules, 
single  photon  absorption  may  happen  along  singlet-singlet  transitions  from  ground  and/or  lowest 
excited  state  (Fig.  3(B0),(Bi)),  and/or  along  triplet-triplet  transition  (Fig.  3(B2))  from  the  lowest 
triplet  excited  state  which  is  populated  by  electrons  relaxed  along  the  intersystem  crossing  link 
(here  we  do  not  take  into  account  ultrafast  relaxations  happening  from  the  vibronic  intermediate 
states).  TPA  occurs  from  the  ground  state  by  simultaneous  absorbing  two  photons,  and  it 
promotes  electrons  to  the  lowest  excited  singlet  state  (Fig.  3(B2))  with  two  possible  following 
transitions:  singlet-singlet  (Fig.  3(Bi))  or  singlet-triplet  (Fig.  3(B2)  from  Nl  to  A, ).  Finally, 

three-photon  absorption  (3PA)  involves  promotion  of  ground  state  electrons  to  the  lowest  excited 
singlet  state  by  simultaneous  absorption  of  three  photons  (Fig.  3(B4)). 
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Fig.  3(a).  Set  of  five  basics  absorption  diagrams:  B0-B4.  Electronic  states  are  labeled  N0-N4, 
and  the  cr’s  are  absorption  cross-sections.  The  upward  arrows  correspond  to  photo-excitation 
transitions  and  downward  arrow  correspond  to  intersystem  electron  decay  (for  simplicity,  we 
omit  all  other  possible  singlet-singlet  and  triplet-triplet  states  electron  decays).  Inset  in  the 
bottom  right  comer  depicts  representation  of  a  TPA-ESA  material  within  our  framework  of 
building  basic  blocks  B0-B4. 
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Fig.  3  (b).  Set  of  three  basic  stimulation  emission  diagrams  B5,  B6,  and  B7.  Similar  to  Fig.  3  (a) 
except  that  the  downward  arrow  corresponds  to  reverse  absorption  in  which  one  (or,  a  >  1 ) 
photon(s)  is  (are)  emitted  from  Ni  to  N0:  number  of  photons  reversely  absorbed  is  specified  as  a 
superscript  of  the  absorption  cross-section  with  negative  value  -a  (this  distinction  is  made  due 
to  the  fact  that  in  many  cases  cross-sections  of  absorption  and  stimulation  emission  are  not  the 
same,  aTPA  *  a^A). 

Diagrams  B0-B7  in  Fig.  3  may  be  considered  as  basic  (computational)  building  blocks  for 
describing  such  materials.  By  combining  “related”  diagrams  one  can  get  description  of  many 
existing  nonlinear  materials.  For  example,  one  can  describe  the  five-level  model  ofC60  -toluene 

solution43  -  nonlinear  RSA  material  -  by  merging  diagrams  B0uBiUB2;  or  describe  the  five- 
level  model  of  a  chromophore  from  AFX  group  "  exhibiting  TPA  assisted  excited  state 
absorption  by  combining  diagrams  B1UB2UB3  which  is  depicted  on  Fig.  3(TPA+ESA).  We 
will  call  nonlinear  material  generic  if  its  energy  diagram  may  be  composed  from  our  set  of  basic 
diagrams  B0-B7. 

A  valid  diagram  composed  from  the  introduced  set  of  transition  diagrams  uniquely 
defines  the  rate  and  propagation  equations.  According  to  the  derivation  described  in  the 
beginning  of  this  section  and  in  Appendix  A  the  rate  equation  in  the  moving  time  frame 
(z,t  =  t'  -kxz\  for  a  material  which  falls  in  our  generic  category  can  be  expressed  in  the  matrix 
form  written  as 


dN  (z,r,t) 
dt 


D 


£>,  +  Z-r^‘'E'V) 

a=l  aflCO0 


N  (z,r,t). 


(7) 


where  N 


y  y  at 

iV0,iVi  iV 


s  electronic  levels, 

D3  =  D3  (  {(7JpA ,  &3PA  }  )  ,  •  •  • , 


is  a  population  density  vector  function  N(z,r,t)  for  a  system  with 

A AeA(K^!?})>  4-4 {{°TPa’°tp1}) , 

®na=DNa({(J[NA PA,vlNNAfpA})  are  NA+\  constant  .S' x.y  sparse 


matrices  of  decay  rates  ktj ,  single  photon  cr. ,  two-photon  <jtpa  ,  three-photon  a3PA ,  and, 

possibly,  Na  -photon  cr[NApA  molar  absorption  cross-sections  correspondingly,  and  I(z,r,t)  is  a 

function  of  photon  flux  density.  The  propagation  equation  for  such  material  is  given  in  the  vector 
form  by 
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(8) 


dl  ( z,r,t ) 
dz 


•/VB 

■2(0, n’Mz,r,t))lP  (z,r,t) , 

P= i 


with  o,  a2=a2{{(7TPA,(TTp}}),  o,=a,({(j3PA,(j(dp}}) 

^Nb=^Nb{{ct[NaWA,ct^nnJpA))  are  NB  constant  (mostly  sparse)  s  dimensional  vectors  composed 
of  certain  elements  of  corresponding  D„  matrices. 

A  better  understanding  of  the  characteristics  of  the  numerical  scheme  for  Eqs.(7)  and  (8)  can 
be  obtained  by  first  making  it  dimensionless  using  the  following  transformations  Tj  =  zjLdf, 

p  —  i  /  R0 ,  r-t/T0,  I(jj,  p,r)  —  1(7],  p,r)  /  /0 ,  Np(j],p,r)  —  Np(7],p,r)/N  ,  L(y  —  n  Rn  j  X ,  /0 
is  the  peak  intensity.  Here  T0,  R0  are  pulse  width  and  beam  radius  described  by  the  incident 
pulse  shape,  and  N  is  the  total  population  electron  density  of  the  material,  which  is  constant  for 
all  times,  N  =  ^  aNg(Tj,p,  r)  .  The  incident  pulse  can  be  described  by  the  general  formula  given 


by  1(7]  =  0 ,p,  r)  =  /0  / (p,r)  ,  e.g.,  f(p,r)  =  exp(-r2)exp(-/?1)  for  the  standard  Gaussian  shape. 
Using  the  transformations  described  above,  Eqs.  (7)  and  (8)  are  given  by 


dN(r],p,  t) 
dr 


=  T 

-  -*■  A 


Dn  + 


na  d  t 

Ta 


L 


ahco() 


Ia(d,P,r) 


N(7,A  r), 


(9) 


mn'p'x) = ■  n(7,  a  n)  ftO,  (io) 

di]  p=x 

respectively.  The  analysis  involves  the  solution  of  Eqs.  (9)  and  (10).  A  description  of  the 
numerical  method  is  given  in  Appendix  A.  Within  the  framework  of  basic  blocks  B0-B7  the 
constants  NA  and  NB  are  restricted  to  be  Aa  <  3  and  NB  <  3  . 


3,1.3  Contributions  of  Each  Electronic  Level  to  the  Total  Absorption 

We  introduce  a  new  feature  into  our  numerical  scheme  that  characterizes  the  materials  from 
our  generic  category,  enabling  an  investigator  to  obtain  a  more  systematic  analysis  of  the 
material’s  nonlinear  properties.  This  feature  is  the  exact  contributions  from  each  electronic  level 
to  the  total  absorption.  Usually  such  contributions  are  roughly  deduced  from  the  dynamics  of  the 
population  densities  of  electronic  levels.  However,  relative  contributions  do  not  always  match 
the  corresponding  relative  population  densities,  so  that  the  approximations  of  such  contributions 
solely  based  on  the  population  densities  may  be  too  far  from  the  actual  values.  Our  numerical 
solution  allows  us  for  any  propagation  distance,  radius,  and  time  step  to  obtain  both  the 
population  density  values  and  intensity  absorption  due  to  each  electronic  level,  so  that  we 
calculate  precisely  the  relative  contributions  to  the  total  absorption.  Appendix  B  gives  the 
detailed  mathematical  analysis  of  the  contribution  of  each  electronic  level  to  the  total  absorption. 

3.2.  Comparison  with  Experimental  Data 

We  selected  a  group  of  nonlinear  materials,  and  validated  our  numerical  scheme  against  their 
transmittance  data  measured  in  different  lasing  conditions.  Our  choice  of  materials  was  mainly 
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motivated  by  a  desire  to  choose  key  representatives  from  the  possible  generic  materials  described 
here.  Among  the  single  photon  absorbers,  we  picked  C60  -toluene  solution,43  which  exhibits 

reverse  saturable  absorption;  for  the  two-photon  absorbers,  we  chose  AFX  chromophore 
AF455; 14-31  ,  and  we  selected  PPAI  dye44  for  the  3PA  materials.  These  materials  attracted  a  lot  of 
attention  in  recent  years  because  of  their  promising  nonlinear  optical  properties. 

Comparison  of  theory  with  experiments  requires  a  definition  of  f(p,r) .  The  analytical 
and  numerical  techniques  explained  above  assumed  a  Gaussian  function  for  f(p,  r) ,  which  well 
describes  the  pulse  shapes  from  the  lasers  used  in  the  experiments.  Therefore  the  incident  laser 
intensity  is  given  by 


7(77  =  0 ,p,  r)  =  70  exp(-  t2)  exp(-/r ) . 


(11) 


3.2.1,  Single  photon  absorbers 

The  first  chosen  nonlinear  material,  C60,  belongs  to  the  group  of  reverse  saturable  absorbers  -  the 
materials  with  ESA  cross  sections  being  much  higher  than  that  of  the  ground  state  -  with  the 
energy  diagram  B0  u  B2  u  B3,  which  uniquely  defines  the  vectors  and  matrices  of  the  coefficients 
for  the  rate/propagation  Eqs.  (7)  and  (8)given  by: 
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(12) 


where  all  the  coefficients  as  well  as  experimental  parameters,  due  McLean  et  al .43,  are  collected 
in  Table  1. 


Table  1.  Parameters  for  multiphoton  absorbing  materials. 


C6o-toluene 

solution3 

AF455  chromophoreb 

PPAI  dyec 

&3PA  (cm3  /  W2) 
<jtpa  ( cm4  /  GW  ) 
o-0i  (cm2) 
a12  (cm2) 

3.1x  10“18  e 

1.6xl0“17 

0.5xlO“2od 

1.68xl0“17 

3.2  x  10“21 

0-34  (cm2) 

1.4xl0“17 

17.1xl0“17 

Ko  (ns) 

32.5 

2.72 

K'  (ps) 

1.0 

1.66  f 

Ki  (ns) 

1.35 

45.3 
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K  (os) 

Kl  (ps) 

40.0  g 

1.0f 

0.368 

10.0  g 

zmax  (mm) 

1.0 

1.0 

10.0 

^df 

0.37306 

2.24 

0.8059 

Nt  ( nm  3 ) 

1.559xl0“3 

h 

0.012  1 

0.596xl0"3  k 

Energy  levels 
diagram 

B0uBiU 

b2 

BiUB2uB3 

b4 

R0  ( jum ) 

16.687  m 

9.25 

53.315 

T0  (ns) 

9.608 

2.71 

2.474x  10“2 

(nm) 

532 

800 

1064 

aMost  of  the  parameters  were  taken  from  Ref.  45,  experimental  parameters  are  described  in  Ref. 
43 

b  Material  and  measurement  parameters  are  described  in  the  experimental  section  of  Ref.  3 1 
c  Ref.  44 
d  Refs.  32,45 
e  Ref.  43 
f  Ref.  14 
g  Ref.  28 

h  Corresponds  to  2.59mM  solution  of  C60  in  toluene 
1  Corresponds  to  0.02M 
k  Corresponds  to  0.99mM  of  the  dye  in  DMSO 

m  All  laser  parameters  are  deduced  from  the  corresponding  original  parameters  by  using  Eq.  (11) 


We  used  the  iterating  scheme  given  by  Eqs.  (A5)  and  (A6)  (See  Appendix  A.)  with 
and  achieved  an  excellent  match  with  the  measured  transmittance  data42  within  the  input 


range 


3.6x10  , 5.0 


J/cm2  ( Fig.  4(a)). 


K  =  2, 
fluence 
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Fig.  4.  Energy  transmittance  TE  as  a  function  of  input  energy.  The  measured  data  is  shown  by 

the  •  (dot)  symbol,  the  solid  line  is  our  numerical  calculation;  the  dashed  line  is  the  original 
solution,  (a)  Comparing  with  the  numerical  solution  of  McLean  et  a/.43  for  C6o,  ns;  (b) 

n  i 

comparing  with  the  analytical  calculation  of  Sutherland  et  al.  AF455,  ns;  (c)  comparing  with 
the  analytical  calculation  of  Wang  et  al.44  PPAI,  ps.  All  the  parameters  are  given  in  the  Table  1. 


Because  the  analytical  solution  of  Eqs.  (7)  and  (8)  in  Ref.  42  is  less  accurate  in  describing 
the  evolution  of  the  population  densities  at  high  fluence  inputs  than  our  complete  numerical 
integration,  our  solution  is  noticeably  better  for  fluence  values  above  1.4  J/cm2 .  Neither  solution, 
however,  matches  the  measured  data  above  5.0  J/cm"  (McLean  et.  al.  noticed  experimentally 
damage  at  high  input  energy).  According  to  the  values  of  the  individual  contributions  collected 
in  Table  2  (see  definitions  in  Appendix  B),  the  nonlinear  transmittance  of  C6o  in  the  ns  regime  is 
mostly  attributed  to  the  lowest  triplet-triplet  state  absorption  varying  from  71%  to  99.8% 
depending  on  the  input  energy  value. 
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Table  2.  Individual  contributions  to  the  total  absorption  of  nanosecond  pulses  to  the  electronic 
states  of  C6o- 


Levels) 

®in(J/cm2) 

0.51  a 

S  4.5  pj) 

(* 

2.05 

17.9  pj) 

(« 

14.1 

123.0  pj) 

(hr 

79.2 

23.1 

10.4 

50.3 

4.0 

0.3 

15.9 

0.0 

0.0 

(A) 

3.8 

5.1 

0.7 

8.5 

3.6 

0.1 

9.8 

0.1 

0.0 

(a) 

16.9 

71.7 

88.8 

41.0 

92.3 

99.5 

74.2 

99.8 

99.9 

(<?)c 

0.06 

0.13 

0.17 

0.08 

0.19 

0.21 

0.14 

0.21 

0.21 

a  For  a  given  input  fluence  value,  we  calculated  three  values  (left,  middle,  and  right  columns) 
corresponding  to  the  relative  contributions  to  the  absorption  averaged  within  every  element  of 
the  T  -range  partition:  [-1,-1/ 3)  ,[-1/3, 1/3) ,  and  [1/3,1]. 
b  Relative  contribution  to  the  absorption  for  level  0 
c  Averaged  intensity  decay  values. 


But  in  addition  to  N3,  the  ground  level  contributes  -23%  to  the  total  absorption  at  low  input 

intensities,  and  this  contribution  is  lost  at  higher  intensities  because  of  the  fast  bleaching  of  the 
ground  level  (see  dynamics  of  the  population  densities  and  individual  contributions  in  Fig.  5). 
This  explains  the  excellent  match  between  calculations  and  experiments  below  1.4  J/cm2 , 
acceptable  match  below  5.0  J/cm2 ,  but  poor  match  of  both  methods  above  5.0  J/cm2.  Material 
degradation  at  high  input  intensities  is  another  possible  reason  for  the  divergence  of  numerical 
solutions,  based  on  integration  of  the  rate -propagation  equations,  from  the  measured  data. 
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Incident  Intensity 


Fig  5.  C6o  nanosecond  regime.  Evolution  of  population  densities  (left  column)  and  absolute 
contributions  (right  column)  due  different  active  electronic  levels  superimposed  with  the  total 
absorption  at  the  entrance  of  the  slab,  lj  =  0 ,  and  at  p  =  0  for  different  incident  fluence  values: 

chin  =  {0.51,  2.05, 14.1}  given  in  [j/cm2  J  (from  top  to  bottom,  correspondingly). 

3.2.2.  Two-photon  absorbers  with  excited  state  absorption 

A  theoretical  basis  for  and  experimental  results  on  nonlinear  transmittance  in  the  ns  regime  of 
D-;r-A  chromophore  AF455  -  a  material  which  exhibits  two-photon  assisted  ESA  -  was 

T  1 

presented  by  Sutherland  et  al.  Within  our  classification,  its  energy  levels  diagram  is  composed 
of  TPA  and  single  photon  absorption  basic  blocks:  B1UB2U  B3,  and  parameters  given  by: 

o,  =[0,crs,a-7.,0,0],  o2  =  [oy^, 0,0, 0,0],  <r3  =  0 
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0 
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l  0 

0 

0 

0 

0; 

where  the  same  matrix  D0 

as  for  C 

60  is 

used.  We  ran 

our 

technique  against  the  experimental 

3 1 

parameters  of  Sutherland  et  al.  (summarized  in  Table  1),  found  excellent  agreement  with  the 
measured  data,  and  observed  nearly  total  coincidence  with  the  previous  analytical  solution  as 
shown  in  Fig.  4(b). 

According  to  the  plots  of  the  population  densities  shown  in  Figs.  6(a),  (b)  for  two  chosen 
input  energy  values,  the  depletion  of  the  ground  level  is  almost  negligible.  It  confirms  the 
validity  of  their  steady-state  approximation  for  the  population  densities  in  ns  regime  that 
explains  the  excellent  match  of  our  numerical  solution  with  their  analytical  solution.  However, 
for  higher  input  energies  this  assumption  fails,  because  there  is  a  considerable  increase  of  the 
population  densities  N{  and  /V,  as  shown  in  Fig.  6(c)  for  the  ns  regime,  whereas  in  fs  regime  the 

ground  state  is  depleted  rapidly,  as  shown  in  Fig.  9(a).  Therefore,  our  numerical  solution  is  more 
accurate  for  higher  input  intensities  as  it  encompasses  the  dynamics  of  population  densities. 
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Incident  Intensity 


Fig.  6.  AF455  nanosecond  regime.  Evolution  of  population  densities  (left  column)  and  absolute 
contributions  (right  column)  due  different  active  electronic  levels  superimposed  with  the  total 
absorption  at  the  entrance  of  the  slab,  rj  =  0 ,  and  at  p  =  0  for  different  incident  energy  values: 

Em  =  jl7pJ,  93pJ,  0.33mJ}  (from  top  to  bottom,  correspondingly). 


An  analysis  of  the  absolute  contributions  from  different  electronic  levels,  whose  plots  are 
depicted  in  Figs.  6(d)-(f),  helps  one  to  better  understand  the  absorption  within  the  AF455  sample 
for  ns  pulses.  The  major  absorption  occurs  in  the  central-to-right  part  of  the  propagating  pulse, 
and  results  in  advancing  the  leading  edge  of  the  pulse  as  it  propagations  through  the  sample.  For 
relatively  small  input  energies,  the  ground  level  is  the  main  contributor  to  the  total  absorption 
with  -54%  of  average  contribution  inside  the  central  part  of  the  pulse  (for  details  see  Table  2 
which  presents  the  individual  contributions  averaged  over  the  leading,  central,  and  trailing  parts 
of  the  propagating  pulse  at  t)  =  0  for  three  input  energy  values).  However,  the  contribution  of  the 
ground  state  TPA  for  higher  energies  (Table  3,  columns  2  and  3)  falls  abruptly,  and  the  main 
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contributors  are  ESA  of  the  lowest  singlet-singlet  and  the  lowest  triplet-triplet  states  with 
corresponding  individual  contributions  of  -40%  and  ~40.9%  for  the  central  part  of  the  pulse, 
respectively. 

Table  3.  Individual  contributions  to  the  total  absorption  of  nanosecond  pulses  to  the  electronic 
states  of  AF455. 


Fevels\ 

£„>J) 

17.0 

93.0 

330.0 

( Po > 

92.8 

54.1 

10.2 

70.5 

18.8 

2.1 

40.9 

7.3 

0.7 

(A) 

5.4 

21.2 

6.3 

22.7 

40.2 

7.0 

47.0 

51.3 

7.5 

(a) 

1.7 

24.5 

83.4 

6.7 

40.9 

90.8 

11.9 

41.3 

91.7 

(<?} 

1.4e- 

1.5e- 

1.3e- 

l.Oe- 

6.0e- 

4 

3 

3 

3 

0.02 

0.03 

3 

0.16 

0.26 

3.2.3.  Three-photon  absorbers 

With  the  advent  of  shorter  pulsed  and  higher  intensity  lasers,  three  or  more  photon  processes  can 
now  be  measured.44  These  materials  may  be  of  interest  in  areas  such  as  laser  micro-or  nano¬ 
machining  because  of  their  ability  to  reduce  the  feature  size  of  structures  in  substrate  materials. 
In  contrast  to  the  previous  nonlinear  materials,  there  is  little  investigation  of  these  materials  such 
as  the  recently  discovered  3PA  dye,  abbreviated  as  PPAI.44  The  initial  experimental  investigation 
of  this  material  shows  a  high  three-photon  absorption  cross-section  in  the  near  infrared  region,  so 
that  all  the  parametric  vectors  and  matrices  will  be  zero,  excluding  the  following  two: 


<*3  =[<t3pa>0]>  A 


PA  0 
V  ^3  PA  0, 


(14) 


and  the  energy  level  diagram  consists  of  only  one  building  block  B4  (the  setup  parameters  are  in 
Table  1).  Our  numerical  solution  matches  extremely  well  with  the  transmittance  measurements 
published  in  Ref.  44  (see  transmittance  curve  in  Fig.  4(c)).  Moreover,  it  is  even  closer  to  the 
measurements  than  the  previously  published  analytic  solution  of  Wang  et.  al.44  in  the  high  input 
intensity  region.  One  obvious  reason  is  that  their  analytical  model  assumes  constant  population 
density  of  the  ground  level.  This  assumption  is  reasonable  in  the  low  intensity  region,  which  is 
demonstrated  by  Fig. 7(a).  But,  in  the  high  intensity  region,  according  to  our  time-resolved 
population  density  solution  depicted  on  Fig.  7(b),  this  assumption  breaks  down  (N0  drops  more 
than  a  half  in  the  second  part  of  the  pulse),  so  that  our  solution  is  more  complete  and  reliable. 
Dynamics  of  the  population  density  /V,  in  Fig.  7(b)  also  confirms  their  supposition  that  excited 
state  absorption  may  significantly  affect  the  transmittance  at  high  input  intensities. 
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Fig.  7.  PPAI  picosecond  regime.  Evolution  of  the  population  densities  at  the  entrance  of  the 
slab,  rj  =  0,  and  at  p  =  0  for  two  incident  intensity  values:  7in  =  {16.9, 204.5}  given  in 

j^GW/cm2  J  ,  which  correspond  to  Em  =  {66.4  pj,  0.8  mj}  . 

3.3.  Femtosecond  Duration  Pulse  Transmission- AF455 
3.3.1.  Predictions  for  AF455 

Unlike  many  previous  numerical  or  analytical  techniques,  our  detailed  analysis  enables  us  to 
investigate  many  cases  of  nonlinear  absorption  including  situations  that  have  not  been  measured, 
particularly  for  short-pulsed  lasers.  The  excellent  agreement  of  our  calculations  with 
experimental  results  gives  us  confidence  to  use  it  as  a  predictive  tool  for  experiments  which  have 
not  yet  been  perfonned.  Since  the  investigation  of  ultrashort  pulses  with  nonlinear  materials  is 
of  increasing  importance  and  the  short  pulse  behavior  of  AF455  has  not  been  quantified  yet,  we 
applied  our  method  using  a  fs  laser  pulse  to  predict  its  effects  on  the  nonlinear  transmission  of 

AF455.  We  assumed  a  Gaussian  pulse  withf?0  =5.0  pm,  T0  =  144.0  fs,  and  \o  =  800 nm  and 
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obtained  -0.002  transmittance  at  the  exit  of  the  AF455  slab  ( zmax  =  3.0  Ldf  ~  0.412  mm)  for  input 
energies  — 10.0  //J  (See  transmittance  plot  for  different  input  energy  values  in  Fig.  8.).  This 
transmittance  value  is  almost  two  orders  of  magnitude  smaller  than  that  obtained  by  using  ns 
pulses.  By  increasing  the  slab  thickness,  one  can  hypothetically  achieve  even  less  transmittance, 
which  makes  AF455  an  ideal  candidate  for  high  intensity  nonlinear  absorbers.  With  further 
validation,  our  numerical  method  may  enable  one  to  substitute  calculations  for  many  laboratory 
experiments,  which  are  both  costly  and  time-consuming. 


Fig.  8.  Transmission  as  a  function  of  input  energy  for  AF455  in  a  0.41  mm  slab  for  femtosecond 
pulses;  (solid  line)  using  integrated  value  (dashed  line)  using  peak  values. 


As  stated  previously  our  technique  also  enables  us  to  examine  the  population  dynamics 
and  the  contribution  of  each  electronic  level  to  the  absorption.  Fig.  9(a)  illustrates  the  evolution 
of  electron  population  densities  at  the  sample  entrance  for  an  incident  energy  of  6.6  //J  .  Because 
of  the  high  energy  of  the  pulse,  the  ground  level  is  quickly  depleted  and  does  not  repopulate  due 
to  slow  decay  rates  kx  0and  k3  0 .  We  calculated  its  individual  contributions  to  the  total  absorption, 

and  their  averaged  values  are  -95%  within  the  first  part  of  pulse,  -48%  in  the  second  part,  and 
falls  to  less  than  10‘2%  at  the  final  part.  The  depopulation  of  the  ground  level  increases  the 
electron  density  of  the  singlet-excited  state  A,  followed  by  fast  excitation  to  N-, ,  so  that  in  the 
second  half  of  the  pulse  its  individual  contribution  is  -51%.  However,  as  can  be  seen  from  the 
relative  magnitudes  of  q0  and  qx  in  Fig.  9(b),  the  TPA  is  more  important  than  ESA  (to  be  more 
precise,  the  total  absorption  g  ,  defined  by  Eq.  (B3)  (See  Appendix  B),  averaged  in  the  first  part 
is  -0.22,  whereas  in  the  second  part  is  -0.08),  which  makes  AF455  mainly  a  TPA  material  in  the 
fs  regime  as  oppose  to  its  TPA  assisted  ESA  behavior  in  ns  regime  (refer  to  Table  3).  Since  the 
intersystem  crossing  time  1  /  k]  3  is  on  the  order  of  nanoseconds,  the  lowest  triplet-triplet  state  A, 

does  not  significantly  contribute  to  the  total  absorption:  its  individual  contribution  is  less  than  10' 
4%  [therefore,  we  do  not  show  it  in  Fig.  9(b)].  From  our  numerical  calculation,  we  find  that  the 
maximum  values  of  N3  and  N4  are  only  about  10  ”6 . 
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Fig.  9.  AF455  femtosecond  regime,  (a)  Evolution  of  population  densities  at  the  entrance  of  the 
0.412  mm  slab  for  144.0  fs  pulse  at  p  =  0  ;  (b)  contributions  of  active  electron  levels  to  the 
absorption  superimposed  with  total  intensity  absorption. 

3.3.2.  Comparison  with  experiments  on  AF455 

We  demonstrated  how  a  conventional  five-level  model  of  AF455  may  be  modified  within  our 
framework  to  a  new  six-level  model,  which  enables  better  agreement  with  transmittance 
measurements  at  high  input  energies.  Using  experimental  femtosecond  transmission  data  on 
AF455  obtained  from  the  Air  Force  Research  Laboratory  at  Wright-Patterson  AFB,  we 
calculated  transmittance  of  laser  light  with  a  wavelength  of  780  nm  and  a  pulse  width  of  175  fs 
through  a  sample  of  AF455,  with  molecular  concentration  A0=1.85xl0ls  cm'3,  by  solving  the 

equations  for  the  five-level  model.  Figure  10  shows  a  good  fit  of  the  five-level  model  to  the 
measurements  in  the  range  from  about  0.01  pj  to  about  2  pj  using  the  following  parameters  '  : 

<r12=  1.68x1  O' 17  cm2,  <t34=  17.1xl0'17cm2,  <JTPA  =  0.51xl0'20  cm4/GW,  k101=2.72ns,  k^  =  45.3ns, 
&2i1=l-66ps,  U,,1  =368ns,  k431  =  10ns.  However,  we  can  see  from  Fig.  10  (insert)  that  at  energies 
exceeding  5  pj  the  calculations  diverge  from  the  experimental  data. 
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Fig.  10.  Experimental  energy  transmittance  of  TPA  chromophore  AF-455  (points)  and  calculated 
results  for  a  new  6-level  model  (solid  curve)  versus  5 -level  model  (dashed  curve).  175fs  (sech2 
FWHM),  780nm,  plane-polarized  pulses  from  a  TiS  regenerative  amplifier  (Clark-MXR  CPA- 
2010)  were  focused  with  either  a  750mm  or  a  250mm  lens  into  a  5.1mm  thick  PMMA  slab 
doped  with  TPA  chromophore  AF-455  at  a  concentration  of  1.85x10  cm'  .  The  front  face  of 
the  sample  was  positioned  at  the  focal  plane  of  the  lens.  The  energy  of  the  pulses  was  varied 
using  a  variable  attenuator  consisting  of  a  rotatable  waveplate  and  fixed  polarizer.  The  range  of 
input  energies  in  this  configuration  was  0.04  to  38.9  pj  .  Pulse  energies  were  measured 
immediately  before  and  after  the  sample  and  the  ratio  of  transmitted  to  incident  pulse  energy 
computed.  The  inset  image  demonstrates  the  superiority  of  the  6-level  model  over  the  existing  5- 
level  model  within  a  range  of  high  input  intensities.  Data  are  courtesy  of  the  Air  Force  Research 
Laboratory. 


We  suggested  adding  a  new  sixth  level  as  shown  in  Fig.  1 1  (5-level  +  6-level).  Then 
using  the  parameters  =  160  fs  and  <t25=0.78x10‘17  cm2,  we  achieved  excellent  agreement 
between  the  calculation  and  experiment. 
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(a)  (b) 


Fig.  11.  Extending  an  existing  (5-level)  model  of  AF455  to  a  new  (6-level)  model  (a)  Five-level 
Jablonski  diagram,  (b)  Basic  Block  diagram  extending  an  existing  (5-level)  model  of  AF455  to  a 
new  (6-level)  model  by  introducing  a  new  excited  energy  level  within  the  singlet  state  group 

Once  a  new  hypothetical  level  has  been  introduced  and  tested  by  our  technique,  one  can 
try  to  validate  it  by  other  techniques,  such  as  z-scan  or  pump-probe  measurements.  We  believe 
that  adding  new  energy  levels  will  be  more  and  more  important  in  the  future  when  higher  energy 
lasers  become  available. 

3.4  General  Numerical  Method  (Electric  Field) 

3.4.1.  Derivation  of  the  theory 

It  is  convenient  to  express  the  electric  field  in  terms  of  a  normalized  function: 

E(z,  r,  t )  =  Q(z,  r,  t)Q'0 ,  Q]  =  s0Q'02 .  (15) 

Using  the  transformations  provided  above  the  electric  field  can  be  written  in  dimensionless 
parameters: 

E(ti,p,t)  =  Q(ti,p,t)Qo.  (16) 

To  include  diffraction  effects  in  the  model,  the  intensity  is  expressed  in  terms  of  a  classical 
complex  electric  field,  Ec ,  provided  in  Eq.  (2),  as  follows: 

I=e„,w„{E2c)  (17) 

The  system  of  equations  can  be  made  dimensionless  using  the  transformations  provided  here 
together  with  the  following  relationship: 
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(18) 


(• E„E,)  =  2a2(Q,Q)  =  -Ql(Q,Q ), 

£0 

which  leads  to  the  following  expression  for  the  field  intensity: 

/(/7,  p,  r)  =  2 nc0Ql  (Q(t),  p,  r),  Q{ij,  p,  r))  =  2nc0Ql  Q(ij,  p,  r) , 

with 


Q-(Q,Q). 

Using  the  transformations  described  above,  the  rate  equation  becomes: 


dN(rj,p,r) 

dr 


=  T 


A  +T,-^-(2nc0Q;)°  Q‘(<l,p,r) 


y=]  CCtlCQf 


N  (r/,p,T). 


Eq.  (21)  can  be  expressed  in  a  more  compact  form  by  using  the  following  substitutions: 


D  =  T  D  D  =T 

M)  —  10Ly0  ’  a  —  10 

which  leads  to  the  following  rate  equation: 

dN(rj,p,r) 


D„ 


aft  (Or 


pnc„' aJ)“  ■ 


dr 


D0  +  T.Dr.Q°(t.P.  0 


a= 1 


N  (tj,P,t), 


(19) 

(20) 

(21) 

(22) 

(23) 


3.4.2  Description  of  the  rate  and  propagation  operators 

A  rate  operator  T  is  introduced  as 


T(^r)  =  T,  +  £Tf(7,r),  (24) 

a= 1 

with  Y0  =  Z)0  and  =  DaQa  {ij,p,v) .  This  leads  to  the  following  compact  form  of  the 

rate  equation: 


dN(rj,p,z) 

dr 


=  T(7,r)oN(7,yO,r), 


(25) 


Using  the  transformation  provided  in  Eq.  (19),  a  dimensionless  propagation  equation,  which  does 
not  contain  a  tenn  to  account  for  diffraction  effects  can  be  written  as: 


mn,P,0\LN 

drj  [ 

or,  equivalently,  as: 


■N(//, p,T))(2nc0Ql)P  '  Qp~\tj,p,t ) 


,/M 


-cLdf\Q(ri,p,T),  (26) 


^^drj’ 1 ^  1  (J1,  A r) ~ C i r)  ’ 


(27) 


p= i 


25 


by  using  the  substitutions  g  p  =  LdfN (lnc0Ql  ^  g  p  and  c  =  cLdf  .  Eq.  (27)  can  be  rewritten  in  an 

operator  form  by  introducing  absorption  operators  ,  O  ) ,  and  Oi  to  account  for  absorption 
due  to  the  intensity,  electronic  density,  and  the  linear  absorptions,  respectively,  where 
0^(77,  r)  =  Gp  ■  N(77,  p,  r) ,  0^(77,  t)  =  Q^~\rj,p,  r)  ,  and  O  L=c.  The  resulting  equations  can  be 
expressed  as: 


® (7. T)  =  “X ‘ N(7> A 0) 2^' (7,  =  (7,  r)0 J (7, r) - , 


>3=1 


7?=1 


and 


dQ(j),p,  r) 
drj 


(28) 


(29) 


An  absorption  operator  that  accounts  for  diffraction  effects  (e.g.,  a  diffraction  operator)  can  be 
written  as: 


^(P7)  =  t 


1-/7  — 


p’ 


with 


h  = 


cOqTq 


A  Fourier  transfonn  of  this  operator  has  the  following  form: 


Vf(p,co)^(\-bcopl, 


F 


8n 

Kdo)n  j 


=  (-i  0))\ 


(30) 


(31) 


(32) 


By  adding  together  the  nonlinear  and  linear  absorption  terms,  the  propagation  equation  can  be 
expressed  in  tenns  of  the  electric  field  as: 

dQ{n’p’T)  =  {®(7. 0  +  'Tip, r)( « Q(Jh P,  r) .  (33) 

dr/ 

This  propagation  equation  includes  complex  values  of  Q(rj,p,z) ,  whereas  the  rate  equation 
provided  in  Eq.  (25)  can  be  expressed  in  terms  of  real  vectors  N(rj,p,r)  and  real  values  of 
Q(r/,p,r). 


3.5  Relationship  between  the  Building  Blocks/decay  rates  and  the  Coupled 
Propagation/Rate  Equations 

A  detailed  description  for  a  more  general  set  of  building  blocks  is  provided  by  describing 
associated  schematic  energy  diagrams,  descriptive  notations,  and  a  technique  for  building  a 
coupled  system  of  rate-propagation  equations  based  on  an  energy  diagram  descriptor,  Jmaterial , 
which  is  expressed  in  terms  of  individual  building  blocks..  Two  types  of  computational  building 
blocks  is  used  to  describe  generic  materials:  absorption  blocks  and  inverse  absorption  blocks. 
The  absorption  blocks  are  associated  with  transitions  from  a  lower  energy  level  to  a  higher  one, 
whereas  the  inverse  absorption  blocks  are  associated  with  transitions  from  a  higher  energy  level 
to  a  lower  one.  Such  building  blocks  are  used  to  describe  energy  level  diagrams  associated  with 
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a  variety  of  photoactivated  materials.  They  can  also  be  used  to  uniquely  define  terms  in  a 
corresponding  coupled  system  of  rate-propagation  equations  which  can  describe  interactions 
between,  e.g.,  electromagnetic  pulses  and  absorbing  materials. 

3.5.1,  Energy  Level  Diagrams  and  Computational  Building  Blocks 

An  energy  level  diagram  can  be  used  to  develop  a  set  of  computational  building  blocks. 
Energy  levels  can  be  associated  with  electronic  states,  vibrational  and/or  rotational  electronic 
substates  in  atoms,  molecules  and  solids,  and/or  bands  and  discrete  levels  in  semiconductors  and 
metals.  An  energy  level  diagram  can  also  define  possible  electron  and/or  exciton  transitions 
between  such  energy  levels  by  describing  transitions  and  corresponding  absorption  cross- 
sections  or  decay  time  values.  An  energy  level  diagram  may  be  constructed  using  quantum 
calculations,  e.g.,  calculations  based  on  Schrodinger  equations,  experimental  measurements, 
and/or  hypothetical  fonnulations  constructing  a  proposed  diagram  or  adding  elements  to  an 
existing  diagram. 

The  Schrodinger  equation  can  include  all  possible  energy  levels  and/or  states  that  a  material 
system  can  be  in  at  a  particular  moment  in  time  and  a  particular  place  in  space,  and  can  describe 
the  evolution  of  such  material  system.  For  example,  a  particle  of  mass  M  (e.g.,  an  electron) 
having  a  potential  energy  Vp(r,t )  can  be  described  by  a  wavefunction  'E SE(r,t)  that  satisfies 
the  following  Schrodinger  equation: 

Pi2  pvp  (T. 

~  —  (T,  t )  +  VP  (7,  t)^SE  (r ,  t)  =  ih  J .  (34) 

2  M  at 

Systems  containing  many  particles,  e.g.,  atoms,  molecules,  liquids  or  solids,  can  be  described 
using  more  complex  but  similar  equations  as  described  in  Ref.  46.  The  probability  of  locating  a 
particle  within  a  volume  du  about  the  position  r  in  a  time  interval  between  /  and  t  +  dt  can  be 
written  as: 


P(r,t)dvdt  =  \'¥SE(r,t)\  dvdt .  (35) 

Allowable  energy  levels  En  of  the  particle  can  be  obtained  by  neglecting  time-varying 
interactions,  using  the  following  form  of  a  wavefunction,  expressed  as  a  product  of  spatial  and 
temporal  terms: 


'PffiCAt)  =  ^(r)exp[i(En /h)t] ,  (36) 

The  spatial  term  t//(f )  can  be  described  by  the  following  time-independent  Schrodinger 
equation: 

~  VV(0+ r,(rM?)  =  £>(?)  ■  (37) 

2  M 

and  the  energy  levels  E"  can  be  detennined  as  eigenvalues  of  the  operator  equation 

Hyf{r)  =  Enyf{r),  (38) 
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where  the  operator  H  can  be  obtained  from  Eq.  (37)  as 

*  ft2 

H  = - V2  +  Vf(F).  (39) 

2  MP 

The  size  of  the  energy  level  diagram  containing  most  or  all  possible  energy  levels  may  be 
infinite,  so  that  a  reasonable  simplification  should  be  used  to  limit  the  number  of  levels  contained 
there.  For  example,  the  size  of  an  energy  level  diagram  can  be  reduced  by  including  only  such 
energy  levels  which  may  be  important  in  describing  the  evolution  of  a  particular  system  (or 
which  may  contribute  the  most  to  a  particular  physical  phenomenon  of  interest),  and/or  those 
which  may  be  revealed  or  explained  by  experimental  measurements.  Conventional  energy  level 
diagrams  may  be  constructed  in  this  way,  and  can  provide  an  appropriate  agreement  between 
mathematical  models  based  on  such  diagrams  and  experimental  observations. 

Certain  materials  can  exhibit  measured  characteristics  and/or  behavior  (e.g.,  optical 
transmittance)  which  may  not  fully  agree  with  the  existing  model  that  can  be  limited  to  a 
particular  set  of  energy  levels.  For  example,  additional  energy  levels  may  need  to  be  added  to  an 
existing  model  as  new  experimental  techniques  and/or  conditions  are  provided  such  as  a  new 
laser  frequency,  a  higher  laser  intensity  obtained  using  a  reduced  temporal  pulse  width,  etc. 
Such  further  hypothetical  energy  level  diagrams  can  be  provided  based  on  experience,  judgment 
and/or  comparisons  with  energy  diagrams  of  similar  materials.  A  mathematical  model  based  on 
a  new  energy  level  diagram  can  be  tested  against  the  experimental  measurements.  An  example 
of  this  approach  is  C6o,  which  was  originally  described  using  an  energy  diagram  having  only 
low-lying  singlet  states  to  explain  experimental  absorption  data.  However,  later  experiments  did 
not  agree  with  this  simple  model,  and  an  additional  set  of  triplet  states  was  added  to  the  energy 
level  diagram.  Determinations  based  on  this  additional  energy  level  diagram  exhibited  a  better 
agreement  with  the  experimental  measurements. 

The  energy  levels  can  be  obtained  as  eigenvalues  of  the  Schrodinger  time-independent 
operator  equation  provided  in  Eq.  (38).  There  can  be  an  infinite  number  of  discrete  energy  levels 
(e.g.,  electronic  levels  of  an  isolated  atom)  associated  with  a  material  based  on  a  solution  to  this 
equation.  The  materials  that  include  many  atoms  (e.  g.,  molecules  or  crystals)  can  have 
additional  energy  levels  associated  with  a  certain  electronic  energy  level.  To  simplify 
construction  of  an  energy  level  diagram,  a  set  of  energy  levels  provided  by  the  Schrodinger 
equation  can  be  restricted  to  those  levels  required  to  generate  reasonable  agreement  between 
calculations  and  experimental  measurements. 

3.5.2  Energy  Level  Diagrams 

For  example,  Figs.  12  (a)-(c)  show  energy  levels  which  can  be  used  in  energy  level  diagrams. 
A  certain  vibrational  and/or  rotational  energy  levels  associated  with  it  which  may  be  degenerate 
with  respect  to  an  energy  level  as  shown  in  Fig.  12(b).  The  electronic  energy  level  can  have 
vibrational  levels  associated  therewith,  and  rotational  levels  can  degenerate  with  respect  to  these 
vibrational  levels,  as  shown  in  Fig.  12(c).  Details  of  the  vibrational  energy  levels  v'  and 
rotational  energy  levels  u'.  are  shown  in  Fig.  12(c).  There  may  be  a  small  difference  between  the 

vibrational  and  rotational  energy  levels  as  compared  to  an  energy  difference  between  various 
electronic  levels.  The  vibrational  and  rotational  energy  levels  can  be  referred  to  as  “sub-levels” 
or  as  a  “manifold  of  states”  associated  with  a  particular  electronic  energy  level  or  “state.” 
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Fig.  12.  Energy  level  diagrams,  (a)  single  electronic  level,  (b)  manifold  of  states,  (c)  vibrational 
and  rotational  states. 


A  small,  finite  set  of  energy  levels  may  be  used  to  describe  certain  energy  level  diagrams. 
The  term  “ energy  level”  can  represent  either  an  energy  level  E”  obtained  as  an  eigenvalue  of  the 
operator  equation  provided  in  Eq.(38),  or  a  degenerate  energy  level  as  shown  in  Fig.  12.  An 
energy  level  can  be  included  in  an  energy  level  diagram,  if  such  energy  level  provides  a 
significant  contribution  to  a  macroscopically  observable  physical  phenomenon  such  as  an 
experimental  measurement  or  a  known  absorption  mechanism. 

In  many  systems  that  include  an  absorbing  material  interacting  with  a  laser  operating  in  the 
UV  to  near  IR  spectral  regions,  the  vibrational  and/or  rotational  energy  sublevels  within  a 
manifold  of  states  may  exhibit  very  short  life  times  compared  to  a  temporal  pulse  width  of  the 
laser.  In  such  systems,  electrons  from  higher  energy  substates  can  relax  to  a  lowest  possible 
substate  within  a  manifold  in  a  time  that  is  much  shorter  than  the  duration  of  a  pulse,  so  that  the 
contributions  of  these  higher  energy  substates  to  a  macroscopic  physical  phenomenon  may  be 
ignored. 

3.5.3.  Energy  Level  Criteria 

The  following  criteria  are  used  to  select  and/or  build  energy  levels  associated  with  particular 
materials  to  form  an  energy  level  diagram: 

1.  levels  that  are  reachable  and/or  accessible  by  entities  energetically  promoted  by  one  or 
more  incident  photons; 

2.  levels  that  can  provide  a  significant  contribution  to  a  macroscopic  physical  phenomenon, 
such  as  one  indicated  or  suggested  by  experimental  measurements; 

3.  levels  which  can  be  degenerate  energy  levels  or  combinations  of  degenerate  levels;  and 

4.  levels  provided  based  on  theoretical  predictions. 

According  to  the  first  criterion  above,  energy  levels  may  be  selected,  which  can  be  reached 
and/or  accessed  by  entities  (e.g.,  electrons  or  excitons)  that  may  absorb  incident  photons.  The 
second  criterion  can  allow  a  selection,  from  the  energy  levels  which  may  satisfy  the  first 
criterion,  of  a  smaller  number  of  energy  levels  which,  when  using  a  mathematical  model,  may  be 
used  to  describe  and/or  account  for  macroscopic  phenomena  such  as  measured  physical 
absorption  processes.  The  third  criterion  can  allow  combination  of  certain  energy  levels  to  form 
degenerate  energy  levels  by  a  combination  of  vibrational  and/or  rotational  sublevels  within  a 
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particular  manifold  of  electronic  states.  The  third  criterion  can  also  allow  a  combination  of  such 
degenerate  levels  into  a  further  “joint”  degenerate  energy  level.  The  fourth  criterion  can  allow 
the  building  blocks  to  be  used  to  expand  conventional  energy  level  diagrams  or  to  build  new 
energy  level  diagrams  by  introducing  additional  levels.  This  procedure  can  provide  a  more 
accurate  description  of  absorption  phenomena  in  a  particular  material  without  requiring  a  new 
algorithm  or  numerical  method  to  solve  a  set  of  equations  describing  the  expanded  energy  level 
diagram  for  the  material.  For  example,  when  using  the  lasers  having  greater  intensity  and/or 
spectral  ranges,  additional  terms  may  be  added  to  an  existing  energy  level  diagram  for  a 
particular  material  to  achieve  consistency  between  the  predictions  and  the  experimental 
measurements. 

3.5.4  Absorption/Reverse  Absorption  Block 

An  absorption  block  from  the  set  B0-B5  represents  a  mechanism  that  includes  simultaneous 
absorption  of  a  photons  by  an  entity  (e.g.,  an  electron  or  an  exciton)  and  dislocation  of  an  entity 
(e.g.,  an  electron  or  an  exciton)  from  a  particular  energy  level  N  to  another  energy  level  /V(  . 

We  distinguish  two  types  of  photon  absorption  mechanisms  in  this  framework:,  a  forward 
absorption,  for  which  a  >  0  ,  refers  to  an  absorption  mechanism  in  which  a  photoactivated  entity 
(e.g.,  an  electron  or  an  exciton)  is  promoted  from  a  lower  energy  level  to  a  higher  energy  level, 
and  a  reverse  absorption  which  refers  to  an  absorption  mechanism  in  which  an  electron  relaxs 
from  a  higher  energy  level  to  a  lower  energy  level  by  re-emitting  one  or  more  photons  that  is 
coherent  with  the  incident  light.  This  mechanism,  for  which  a  <  0 ,  describes  the  phenomenon 
of  stimulated  emission. 

For  example,  an  absorption  block  diagram  shown  in  Fig.  13  includes  indices  Sj  ,s2  which  may 
represent  certain  energy  levels  that  an  electron  is  promoted  “from”  and  “to”  respectively  (or 
relaxed  “to”  and  “from”  respectively  if  a  <()).  The  parameter  a  represents  the  number  of 
photons  that  is  absorbed  to  promote  an  electron  from  level  Sj  to  level  s0 ,  where  a  can  be 
negative  if  the  expression  is  used  to  represent  a  reverse  absorption  block  that  includes  a 
relaxation  from  level  s2  to  level  5, .  Optionally,  the  parameter  a  may  be  omitted  if  it  is  equal  to 
one. 


Fig.  13.  Energy  level  diagram  for  absorption  of  a  -photons. 
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3.5.5.  Example  of  absorption  block  descriptor 

The  energy  diagram  for  AF455  can  be  represented  as  the  following  energy  level  diagram 
descriptor,  using  the  computational  building  block  notation  described  here: 


Zap  455  ={KB2>Bi}  (4°) 

where,  B3  represents  a  simultaneous  absorption  of  two  photons  and  promotion  of  one  electron 
from  a  ground  level  to  a  lowest  excited  level  of  the  singlet  state;  B2  represents  an  absorption  of 
one  photon  and  promotion  of  one  electron  from  the  lowest  excited  state  to  the  higher  excited 
state  in  the  triplet  state;  Bi  represents  an  absorption  of  one  photon  and  promotion  of  one  electron 
from  the  first  excited  state  to  the  higher  excited  state  in  the  singlet  state. 

3.5.6.  Stimulated  emission 

An  example  of  stimulated  emission  can  be  described  from  a  two-photon  absorption  (TPA)  state 
with  nonlinear  excited-state  absorption  (ESA).  The  emission  may  be  to  a  conduction  band  or  it 
may  be  a  free  electron  excitation.  A  material  that  can  exhibit  stimulated  emission  effects  and 
which  may  be  described  is  a  Green  Fluorescent  Protein  (GFP)  47 

For  example,  an  electronic  configuration  of  a  material  is  shown  in  the  figure  below  which 
includes  certain  states:  a  ground  state  manifold  No,  a  first  excited  state  manifold  Ni,  which  may 
be  reached  through  a  TPA  event  and  can  be  referred  to  as  a  TPA-state,  and  a  second  excited  state 
N2,  which  can  be  reached  via  a  one  photon  absorption  event  from  the  first  excited  state  manifold. 
The  second  excited  state  can  be  referred  to  as  a  ESA-state,  and  it  may  be  a  conduction  band  or  a 
free  electron  excitation.  A  stimulated  emission  can  occur  between  the  first  excited  state  and  the 
ground  state.  A  spontaneous  emission  can  occurs  between  the  second  excited  state  and  the  first 
excited  state,  or  between  the  second  excited  state  and  the  ground  state. 


(a) 


(b) 
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Fig.  14.  Green  Fluorescent  Protein:  (a)  Jablonski  diagram  of  TPA  and  ESA.  (b)  Building  Block 
diagram  showing  stimulated  emission  and  decay  rates. 

Parameters  kio  and  k7i  shown  in  Fig.  14  can  represent  decay  rates  of  the  TPA-state  and  the 
ESA-state,  respectively,  and  the  parameter  koo  can  represent  a  decay  rate  from  second  excited 
state  1  to  the  ground  state.  The  corresponding  matrices  that  may  be  used  in  the  rate  equation  can 
be  written  as: 
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The  parameter  g10  =  g1  /  g°  may  also  be  used,  where  g1  and  g°  can  represent  an  electron 
degeneracy  of  electronic  levels  1  and  0,  respectively.  In  the  present  example,  all  other  matrices: 
DNa  (ajv^)  =  °;  cTj  =(0,<r£ai,0),  a2  =  (gwcrTPA,aTPA,0) ;  and  all  other  vectors 

(f7[vs]^)  =  0- 

3.5.7.  Absorption  blocks  and  decay  rates 

Each  absorption  building  block  and  corresponding  relaxation  rate  matrix  can  define  a  unique 
corresponding  term  in  the  coupled  system  of  rate-propagation  equations  provided  in  Eqs.  (A5) 
and  (A6),  which  may  be  used  to  describe  propagation  of  an  electric  field  through  a  sample  of  a 
photoactive  material.  The  computational  building  blocks  described  here  are  analogous  or  similar 
to  computational  molecules  that  are  used  to  formulate  multidimensional  finite  difference 
numerical  schemes. 

A  particular  energy  diagram  descriptor  xmat  associated  with  a  photoactivated  material  mat 
can  include  a  component  b  e  %mat .  The  component  b  can  represent  one  of  the  types  of  the 
absorption  building  block  or  the  stimulated  emission  building  block.  Each  computational 
building  block  can  provide  certain  terms  in  a  matrix  Da ,  a  >  0  ,  and  a  vector  a fj  which  appear 
in  the  rate  and  propagation  equations  provided  in  Eqs.  (A5)  and  (A6). 

3.5.7.1  Representation  of  absorption  block/decay  rate  parameters 

For  example,  b  can  represent  an  absorption  building  block,  B0-B7,  which  may  be  associated 
with  an  absorption  parameter  <7**  ,  cr[a]PA ,  or  <y[~“pA,  where  Si  and  S2  are  described  above  and 
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a  >  0 .  The  value  of  this  parameter  can  be  included  at  corresponding  places  in  the  matrix  Da  and 
vector  a g ,  which  can  lead  to  the  following  expressions: 

For  a  =  1  and  as^  ,  the  expressions  are  given  by 


D|[.v|,.V|]  =  D][.sp.S|]-rrv|S2 
D|[.v2,.sl]  =  D|[.v2,.vl]  +  crV2, 

o,bi]  =  (T1[51]  +  cr^2. 


For  a  =  1  and  ,  the  expressions  are  given  by 


and 


-Di[52,52]  D^S2,S2\  <7^ 

Dl[.y,,>s’2]  =  D|[.y|,>s’2]  +  (7s(;'), 


o1[^1]  =  o1[51]-rrs(-1). 


(42) 

(43) 

(44) 


(45) 

(46) 

(47) 


For  a  >  1  and  cr[a]PA ,  when  an  absorption  transition  involves  states  y  and  s2 ,  the  expressions  are 
given  by 


Az[Wi]  =  Da[svsx\-cj[a]PA 
Da\_S2->S\  \  =  Da[s2,sx]  +  <J{a}PA 

G«l>l]  =  a«l>l]  + 


(48) 

(49) 

(50) 


For  a  >  1  and  (j[a“TlA ,  when  an  absorption  transition  involves  states  s2  and  s1 ,  the  expressions  are 
given  by 


and 


Da\_s  2^2]  Da[s  2,S2]  ^[aJPA 
Da  [^1 5  S2  ]  =  Da  IX ;  ^2  ]  +  ^[a]PA 


(51) 

(52) 


(53) 


The  various  decay  rates  for  each  transition  can  placed  in  the  matrix  D0  according  to  Eq.  (9). 


2.5.1. 2.  Examples  of  inclusion  of  parameters  in  matrices  and  vector 

The  matrices  D0,  Da  and  the  vector  o/;  which  appear  in  the  rate  and  propagation  equations 
provided  in  Eqs.  (A5)  and  (A6)  may  be  constructed  for  a  particular  material  using  the  following 
procedure.  For  example,  each  matrix  D0,  Da  and  the  vector  o/;  may  be  defined  initially  to 
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have  all  elements  equal  to  zero.  The  energy  diagram  descriptor  %mat ,  for  a  given  photoactive 
material  can  be  obtained  using  the  procedures  described  here.  Each  computational  building 
block/decay  rate  appearing  in  xmat  can  then  be  used  to  modify  certain  elements  in  the  matrices 

D0 ,  Da  and  the  vector  o/; .  For  example,  for  a  =  1  and  a  ,  Eqs.(42)-(44)  or,  for  a  =  1  and 
<j[~^ ,  Eqs.  (45)-(47)  can  be  applied  for  each  absorption  block  that  appears  in  the  energy  diagram 
descriptor.  As  another  example,  for  a  >  1  and  &[a]PA ,  Eqs.  (48)-(50)  or,  for  a  >  1  and  c^pa , 

Eqs.  (51)-(53)  can  be  applied  for  each  absorption  block  that  appears  in  the  energy  diagram 
descriptor.  This  procedure  can  be  repeated  until  all  component  b  in  the  energy  diagram 
descriptor  (e.g.,  all  absorption  coefficients  and  decay  rates)  have  been  used  to  modify  appropriate 
elements  in  the  matrices  D0 ,  Da  and  the  vector  a fl  using  Eqs.  (42)-(53).  The  resulting  matrices 

and  vector  may  be  provided  in  Eqs.  (A5)  and  (A6),  and  these  equations  may  be  solved  using 
procedures  described  here  to  determine  the  absorption  behavior  of  the  particular  material. 


3.6  Semiconductor  Quantum  Dots 
3.6.1  Fabrication  Techniques 

Semiconductor  quantum  dots  (QDs)  are  often  referred  to  as  artificial  atoms  (See  figure  below), 
which  consist  of  hundreds  to  thousands  of  atoms.  A  few  reviews  of  the  synthesis,  structure  and 
properties  of  nanomaterials  are  given  the  references  48-50.  In  the  1980s  researchers  formulated 
the  concept  of  size  quantization51'53  which  was  later  expanded54’55  As  a  result  of  the  quantum 
confinement,  QDs  exhibit  unique  physical  and  optical  properties  that  are  not  present  in  bulk 
material.  The  terminology  used  for  QDs  evolved  over  several  years  and  is  reviewed  in  ref.  56. 
A  material  is  denoted  as  a  QD  if  the  size  of  the  material  is  smaller  than  the  Bohr  radius  i.e.  the 
size  of  material  is  less  than  the  extension  of  the  electron-hole  pair  in  bulk  material.  In  this  case, 
an  exciton  formed  in  the  QD  has  no  degrees  of  freedom  and  no  motion  is  possible  for  the 
exciton. 


(a)  (b) 


34 


Fig.  15.  Schematic  diagrams  of  a  semiconductor  quantum  dot  (QD):  (a)  QD  alone  where  each 
atom  is  represented  by  the  gold  sphere;  (b)  core/shell  structure  (QD  with  organic  shell  which  is 
represented  by  the  blue  lines). 

3.6.1. 1  Methods  of  colloidal  synthesis 

There  are  common  methods  for  chemical  synthesis  of  QDs:  (1)  Arrested  precipitation  using  a 
stabilizing  agent.  (2)  Reverse  micelles.  An  example  of  the  first  method  is  given  by  the 
formation  of  CdS  QDs  in  which  a  Cd  cupferronate  solution  in  toluene  is  poured  on  the  surface  of 
Na?S  aqueous  solution  to  form  QDs  at  the  interface.  The  QD  size  is  controlled  by  varying  the 
reagent  concentrations,  temperature  and  reaction  duration.  A  common  stabilizer  for  the  reaction 
is  tri-octylphosphine  oxide  (TOPO).  An  example  of  the  second  method  is  given  as  follows: 
CdS  QDs  are  formed  by  preparing  two  mixtures,  water-toluene  and  water-surfactant.  One 
aqueous  phase  contains  the  required  amounts  of  a  Cd  salt;  while,  the  other  one  contains  sodium 
sulfide.  Mixing  these  two  solutions  results  in  micelle  capped  QDs  as  shown  in  part  (b)  of  the 
above  figure. 

3.6.1.2  Core-shell  quantum  dots 

The  surface  to  volume  ratio  in  QDs  plays  a  significant  role.  A  review  of  surface  properties  is 
given  in  ref.  59.  In  order  to  passivate  the  QD  surface,  other  materials  are  used  forming  a  so- 
called  core-shell  QD.  Inorganic  shells  are  often  fabricated  from  another  semiconductor  which 
differs  from  the  core  material.  Common  examples  of  semiconductor-semiconductor  core-shell 
QDs  are  CdS/CdSe  and  CdSe/ZnS.  A  common  synthesis  comprises  a  shell  buildup  by  deposition 
of  cations  and  anions  on  the  surface  of  QDs  by  successive  ion  layer  adsorption  and  reaction60. 
A  semiconductor-semiconductor  core-shell  QD  can  be  thought  of  as  a  solid  solution  of  two 
chalcogenides;  typically,  these  materials  have  higher  efficiency  of  photo  luminescence.  The 
reverse  micelles  also  form  a  core-shell  structure,  but  their  photoluminescence  properties  have  not 
been  studied  in  as  much  detail  as  the  semiconductor-semiconductor  core-shell  QDs. 

3.6.1.3  Quantum  dot  nanostructures 

3.6.1.3.1  Quantum  dots  on  substrates 

In  many  applications  a  planar  structure  is  needed  as  in  photo-electronic  or  electronic  devices. 
Methods  include  molecular  epitaxy  and  metal  organic  chemical  vapor  deposition,  which  are  well 
known  in  semiconductor  laboratories61.  Chemical  methods  can  also  be  used  such  as  self- 
assembled  on  a  surface.  A  monolayer  of  CdS  QDs  on  a  substrate  was  perfonned  by  fabricating 
reverse  micelles  stabilized  with  sodium  bis(2-ethylhexyl)sulfosuccinate  .  Additional  treatment 
with  dithiol  created  a  multilayer  structure  for  the  QDs63.  Layer-by-layer  adsorption  techniques64 
have  also  been  used  to  create  organized  QD  structures. 

3.6.1.3.2  Quantum  dots  in  inorganic  matrices 

The  small  size  of  the  inner  cavities  of  zeolites  allows  for  the  synthesis  of  QDs.  An  example  of 
this  method  for  CdS  and  PbS  was  investigated  in  ref.  65.  In  this  case  the  QD  size  was  restricted 
to  1  mn.  Another  method  involves  mesoporous  silica  gel66  as  a  matrix  for  QD.  An  example  is 
given  by  CdS  QDs  bonded  to  Zn  atoms  of  Zn-doped  material  through  dithiol  bridges  .  These 
structures  yield  a  high  photocatalytic  activity. 
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3.6.1.3.3  Quantum  dots  in  polymers 

QD  in  polymers  is  a  very  active  field.  An  example  is  given  by  modifying  CdS  QDs68  on  the 
surface  of  SiCfi  nanogranules  by  using  a  radical  polymerization  initiator  followed  by 
polymerization  of  the  methyl  methacrylate.  The  polymer  solution  with  the  QDs  was  poured  on 
the  substrate  and  the  solvent  was  removed.  Researchers  believe  that  high  photoluminescence  can 
be  achieved  with  QDs  well  distributed  in  a  polymer  matrix.69 

3.6.2.  Electronic  levels  of  quantum  dots 

A  quantum  dot  (QD)  is  a  quasi-zero-dimensional  object  where  the  carrier  movement  is  restricted 
in  three  dimensions.  The  bulk  crystalline  structure  of  the  semiconductor  is  maintained  in  the 
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QD,  but  they  have  atomic-like  discrete  energy  states  due  to  3D  confinement. 


(a)  (b) 

Fig.  16.  Energy  level  diagrams,  (a)  bulk  semiconductor  (b)  quantum  dot  assuming  single  valence 
bands. 

Carrier  dynamics  in  QDs  are  quite  different  from  bulk  materials  due  to  the  large  energy  level 
spacing  and  enhanced  surface-to-volume  ratio.  State-filling  is  the  main  mechanism  for  optical 
nonlinearities.71"72  Typically  the  effective-mass  model  is  used  to  describe  the  electronic  structure 
of  QDs73.  In  this  description,  electron  states  are  annotated  by  using  a  letter  (/)  to  denote  the 
angular  momentum  of  the  envelop  wave  functions.  Using  spectroscopic  notation,  the  following 
symbols  are  used  S  for  1=0,  P  for  1=1,  D  for  1=2,  etc.  The  degeneracy  of  the  electron  states  are 
2(2 1  +  1)  -fold  degenerate  and  the  lowest  three  electron  states  are  denoted  by  IS,  IP,  ID.  The  hole 
states  are  generally  more  complicated  than  shown  in  the  figure  above  and  are  labeled  similarly  to 
electron  states  with  an  additional  subscript  for  the  total  hole  angular  momentum  ( j  ),  which  have 
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a  degeneracy  of  2  /  + 1  .  The  first  three  hole  states  are \Sm,\Pm, 2S3/2 .  The  interband  selection 
rules  allow  for  transitions  from  nSXh)  hole  states  to  all  S(e)  electron  states,  nPXh )  hole  states 
to  all  P(e)  electron  states,  etc.  where  n  is  the  ordinal  number  of  the  level.  A  schematic  diagram 
is  given  in  the  figure  below. 


Fig.  17.  Energy  level  diagram  and  transitions  for  semiconductor  quantum  dot. 

The  above  figure  shows  five  absorptive  transitions  of  a  single  photon  from  states  IS3/2  to  2S; 
2S3/2  to  2S;  IP3/2  to  IP;  IP1/2  to  IP  and  IS  1/2  to  IP.  The  energies  of  the  excitons  depend  on  the 
relative  energies  of  the  electron  and  hole  states.  The  energies  of  the  exciton  transitions  may  be 
obtained  from  spectroscopic  measurements.  Often  these  spectroscopic  measurements  are 
compared  to  theoretical  calculations  of  the  electron  and  hole  energy  levels.  The  figure  below 
shows  the  theoretical  calculations  for  PbSe  QDs.74  Part  (a)  of  the  figure  shows  the  theoretical 
electron  and  hole  energy  levels  as  a  function  of  QD  size;  while,  part  (b)  shows  the  theoretical 
exciton  energy  states  from  part  (a)  following  transition  rules  and  electron  degeneracy  as 
described  above. 
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Fig.  18.  Electronic  levels  for  PbSe  QDs.  (a)  Theoretical  energy  levels  as  a  function  of  QD  size, 
(b)  Theoretical  transitions  for  four  QD  sizes  with  experimental  transitions  shown  as  dots. 
(Revised  figure  from  Ref.  74.) 

3.6.3.  Generation  and  recombination  of  excitons 

There  are  several  generation/recombination  processes  in  semiconductors76'80.  The  table  below 
describes  the  major  processes.  A  particular  process  known  as  the  Auger  process  is  non-radiative 
and  is  mediated  by  electron-electron  (e-e)  Coulomb  interactions,  which  differ  amongst  atoms, 
bulk  semiconductors  and  QDs.  In  atomic  systems  the  e-e  coupling  is  much  stronger  than  the 
electron-phonon  coupling  which  makes  the  rates  of  Auger  transitions  much  greater  than  the 
radiative  transitions.  Thus  the  decay  of  multi-electron  levels  is  dominated  by  Auger  processes. 
However,  in  bulk  semiconductors  the  efficiency  of  the  Auger  process  is  reduced  because  the  e-e 
interactions  are  smaller.  The  situation  in  QDs  differs  from  the  bulk  in  that  the  quantum 
confinement  enhances  the  e-e  leading  to  increased  Auger  rates.  Additionally,  QDs  have  a  large 
surface  to  volume  ratio  that  leads  to  trap  states  in  the  band  gap  (not  shown  in  the  table).  In  many 
QDs  the  surface  states  are  passivated  by  capping  the  QD. 


Table  4.  Type,  description  and  schematic  diagram  of  various  recombination 
mechanisms  of  excitons76'80. 


Diagram 

Band-to-band:  single  electron 

decay  or  absorption 

Generation  or  recombination: 

Electron  decay  (shown  on  right)  or 
absorption  from  its  state  in  the 
conduction  or  valence  band  and 
occupies  or  creates  the  hole  in  the 
valence  band:  Ephoton  >  Egap 

_ A _  GR 

X 

l 

o 

Trap-assisted 

Electron  and  hole  combine  in  a  trap 
state  in  the  band  gap 

_ A _  r.R 

> 

3 

; - X 

c 
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Auger 

Electron  and  hole  combine  in  a 
band-to-band  transition  and  the 
energy  is  absorbed  by  another 
electron  in  the  conduction  band  or  a 
hole  in  the  valence  band 


CB 


VB 


Impact  ionization 

Generation  mechanism  that  is  the 
counterpart  to  Auger  recombination. 
Occurs  when  an  electron  has  energy 
much  larger  than  the  conduction 
band  edge.  The  excess  energy  is 
generates  an  electron-hole  pair 
through  band-to-band  transition. 


T7 


CB 


VB 


Stimulated  emission 

Population  inversion  exists  and 
incoming  photon  induces  emission 


CB 


VB 


Multiphoton  absorption 

Absorption  of  two  or  more  photons 
simultaneously. 


3 


T 


CB 


VB 


± 


Carrier  multiplication 

More  than  one  exciton  is  created: 
Example  of  bi-exciton  i.e.  creation 
of  two  excitons  on  one  state. 


CO 


CB 


VB 
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Transparency 

For  a  ground  state  containing  two 
electrons,  excitation  of  one  electron 
can  lead  to:  a  photon  excitation  of 
the  ground-state  electron,  which  is 
absorbed.  The  medium  is 

transparent  because  there  is  no  net 
gain  or  loss  of  photons. 

_ 9 _ 9 _  r.R 

X  ' 

c 

^  VB 

Optical  gain 

Optical  gain  may  occur  when  there 
are  more  electrons  in  the  excited 
state  than  in  the  ground  state- 
population  inversion.  If  the  gain 
from  stimulated  emission  is  larger 
than  losses  due  to  absorption  or 
scattering  of  photons,  then  the 
system  may  exhibit  amplified 
stimulated  emission. 

_ 9 _ 9 _  r.R 

X 

3 

uo 

3.6.4.  Three  photoexcitation  regions 

There  are  three  energy  regimes  to  consider:  (a)  E>  Eg,  (b)  E~  Eg,  (c)  E<  Eg  where  E  is  the  energy 
of  the  laser  and  Eg  is  the  energy  of  the  semiconductor  bandgap.  The  table  below  summarizes  the 
various  effects  and  properties,  which  are  then  discussed  below. 


Table  5.  Properties  and  mechanisms  for  three  energy  regimes  in  semiconductor  quantum  dots. 


(a)  E»ES 

(b)  E~ES 

(c)  E<Ea 

Single  photon 

E=Eg 

E<~Eg 

Multi-photon 

Single 

photon 

Multi-photon 

Property: 

one  photon  whose  energy 

exceeds  2  Eg  can  create  multiple 
excitons 

Quantum  Yield:  >>100% 

Property: 
single 
photon 
creates  one 
exciton 

Property: 
multiphotons 
(simultaneously) 
create  one 

exciton 

Property(hypothetical): 
multiple  photons  whose 
total  energy  exceeds  2  Eg 
may  create  multiple 

excitons 

Quantum  Yield:  >>100% 
Experiments  required 

Mechanism:  two  proposed 

Mechanism:  Two-level  model 

Mechanism:  unstudied 

Nozik  et.  al. /4 

Klimov  et.  al.81" 

83 

High  energy 

High  energy 
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exciton 

exciton  decays 

decays  by 

by  sequential 

sequential 

Auger 

phonon 

recombinations: 

emission: 

Noncoherent 

Coherent 

process 

3.6.4.I.  Case  (a)  E»Eg 

In  this  case  a  single  high  energy  photon  can  create  multiple  excitons  resulting  in  high  quantum 
yields  in  excess  of  100%.  The  exact  mechanism  is  not  clear  yet.  However,  there  are  two  main 
groups  involved  in  the  experimental  and  theoretical  analysis.  One  is  at  the  National  Renewable 
Energy  Laboratory/Univ  of  Colorado74  and  the  other  is  at  Los  Alamos  National  Laboratory81.  A 
typical  energy  level  diagram  for  a  semiconductor  quantum  dot  in  the  ground  state  with  electron 
population  is  given  in  the  figure  below.  Both  groups  find  quantum  yields  (QY)  in  excess  of 
1 00%— so  far  up  to  700%.  This  phenomenon  does  not  happen  in  bulk  semiconductors. 


Ground  state 


2  electrons 
6  electrons 


Fig.  19.  Schematic  diagram  of  ground  state  of  a  semiconductor  quantum  dot. 

The  National  Renewable  Energy  Laboratory/Univ  of  Colorado  group  measured  multiple 
exciton  generation  (MEG)  for  single  photon  absorption  in  colloidal  PbSe  and  PbS  QDs.  Using 
transient  absorption  spectroscopy  they  obtained  data  for  intraband  and  interband  probe  energies. 
The  magnitude  of  the  photoinduced  absorption  hange  at  the  band  edge  is  proportional  to  the 
number  of  excitons  created.  The  transients  are  detected  by  probing  either  with  a  band  edge 
probe  pulse  or  a  mid-IR  probe  pulse.  They  introduce  a  model  for  MEG  based  on  the  coherent 
evolution  of  the  quantum  state  from  the  initial  single  exciton  to  a  coherent  superposition  of 
multiple  exciton  states.  The  QD  cooling  via  MEG  occurring  in  PbSe  is  accomplished  by  energy 
loss  through  dephasing  of  the  coherent  multiexciton  state  by  phonon  emission. 

The  Los  Alamos  National  Laboratory  group  also  uses  transient  absorption  spectroscopy  of 
CdSe  QDs.  In  this  case  the  multiple  exciton  formation  is  denoted  as  carrier  multiplication  (CM). 
The  decay  mechanism  is  theorized  to  be  due  to  nonradiative  Auger  recombination.  These 
dynamics  become  increasingly  faster  as  the  number  of  excitons  increases.  Owing  to  strong 
electron-electron  interactions,  a  high-energy  electron  in  the  conduction  band  is  not  stable,  but 
exists  as  a  virtual  electron.  This  virtual  electron  loses  its  energy  through  the  Coulomb  interaction 
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to  a  valence-band  electron.  This  valence-band  electron  is  then  excited  to  the  conduction-band. 
This  process  is  noncoherent.  The  figure  below  shows  the  two  processes. 

MEG  with  phonon  decay  CM  with  Auger  decay  (II) 


Real 

electron 

(black) 


Real 

hole 

(blank) 


rs 


— 

-V 

. \ 
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Fig. 20.  Creation  of  multiple  excitons  with  single  high  energy  photon.  The  green  arrow  represents 
the  high  energy  absorbed  photon,  (a)  Multple  exciton  generation  (MEG)  with  phonon  decay74. 
The  purple  curves  represent  the  phonon  decays,  (b)  Carrier  multiplication  (CM)  with  Auger 
recombination  .  The  blue  line  represents  the  Coulomb  interaction  between  the  virtual  electron 
and  hole.  The  red  dot  represents  a  virtual  electron  and  the  yellow  dots  represent  virtual  holes. 


3.6.4.2.  Case  (b)  E~Eg 

This  is  the  most  commonly  studied  as  it  is  the  easiest.  It  appears  to  follow  a  two-level  model. 


Excite  1  exciton  (low  energy)  Excite  1  exciton  (higher  energy) 


Fig.  2 1 .  Creation  of  a  single  exciton  by  an  absorbed  photon  near  Eg  (a)  Lowest  energy  ( 1  S3/2(h) 
to  lS(e)).  (b)  Higher  energy  (lP3/2(h)  to  lP(e)). 
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In  the  above  figure,  a  single  exciton  is  created  by  a  single  photon.  Depending  on  the  energy  of 
the  incident  photon  different  energy  levels  can  be  excited  and  will  appear  in  the  absorption 
spectrum. 


Excite  1  exciton  (low  energy)  with  2PA 


Excite  1  exciton  (higher  energy)  with  2PA 


(a) 

Fig.  22.  Creation  of  a  single  exciton  by  multiple  (2)  photons  near  Eg  (a)  Lowest  energy  (lS3/2(h) 
to  lS(e)).  (b)  Higher  energy  (lP3/2(h)  to  lP(e)). 

In  the  above  figure,  a  single  exciton  is  created  by  multiple  (e.g.  2)  photons.  Depending  on  the 
energies  of  the  incident  photons  different  energy  levels  can  be  excited  and  the  2PA  can  be 
observed  in  measurements  such  as  z-scan. 


The 
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carrier  absorption,  A,  h  is  the  concentration  of  photoexcited  electron-hole  pairs,  r  is  the 
effective  decay  time,  (3  is  the  2PA  absorption  coefficient. 


3.6.4.3.  Case  (c)  E<  Eg 

As  far  as  we  know  this  case  has  not  been  studied.  It  would  require  a  very  intense  light  source 
and  materials  with  very  large  nonlinear  absorption.  Intriguingly,  if  this  case  behaves  similar  to 
case  (a)  i.e.  single  photon  absorption,  then  when  multi-photon  absorption  exceeds  2Eg  multiple 
excitons  may  be  created.  As  far  as  we  know  this  has  not  been  experimentally  observed,  but  may 
be  present.  This  situation  is  quite  interesting  because  a  very  high  QY  may  be  obtained  in  this 
region. 

4.  Conclusions 

In  this  report,  we  have  described  results  including  a  time-resolved  finite-difference  numerical 
scheme  that  is  more  general  than  traditional  techniques.  The  method  was  validated  against  both 
published  research  and  internal  AFRL  experiments.  Additionally,  we  expanded  our  method  to 
include  electric  field  propagation  and  diffraction  effects.  Finally,  we  investigated  semiconductor 
quantum  dots  in  terms  of  fabrication  techniques,  electronic  energy  levels  and 
absorption/recombination  phenomena. 
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Appendix  A  Description  of  the  Numerical  Method 

We  introduce  a  family  of  identical  2D  grids  Q  indexed  by  radius  p ,  such  that 


n={a.Pj),Pl=jhp),  (Ai) 

«(P,)  =  K(A«,0')),  (A2) 

For  every  p  .  sample,  a  member  from  Q  is  a  pair  of  interleaved  grids  in  the  rj  —  T  parametric 
domain,  such  that  one  grid  is  given  by 

^(j)  =  {(vn+u2,Pj,Ti+m), 

n„+ 1/2  =  (%  +  A/7 / 2)  +  n  A//,  tMI2  =  (r0  +  Ar  / 2)  +  7  Ar} ,  (A3) 

and  is  used  to  sample  the  population  density  N(r/,p,  r) .  Another  grid  given  by 


ni(j)  =  {(TIn>Pj>Ti)>T?n  =r/0  +  n/±ri,Ti  =  r0  +  /Ar}  (A4) 

is  used  to  sample  the  intensity/ (77,  p,r) .  We  integrate  the  resulting  dimensionless  equations 
given  by  Eqs.  (9)  and  (10)  for  a  small  step  sizes  At,  At;  at  Q^(y),  Q7(y)grid  points 

correspondingly  in  two  steps.  In  the  first  step  we  integrate  the  system  of  the  rate  equations  to  find 
a  spatially-resolved  solution  in  the  current  moving  frame  (reference  system  of  the  “pulse  rest”), 

while  keeping  the  average  intensity  —  [ / (rjn , p. , r,)  +  / (/;„  , , y>; , r, ) ]  fixed.  In  the  second  step  we 

solve  the  propagation  equation  for  a  thin  slice  [77, 77  + A/7],  by  using  available  average  population 

density  ^[n(t7„+1/2,/7  .,r,._1/2)  +  N(7„+1/2,/7;,ri+1/2)]  as  an  approximation  for  electronic 

populations  (see  Appendix  B  for  detailed  derivation).  The  resulting  system  of  coupled  equations, 
given  by  Eqs.  (A5)  and  (A6),  still  has  interdependences  between  the  intensity  and  population 
densities;  therefore,  we  do  k  =  \  :K  iterations  to  obtain  a  numerical  solution  of  the  intensity 
function  I[rjn,p,r\  for  a  given  depth  tjn  (for  notational  simplicity  only  the  indices  n,j,  and  i 

will  be  kept  in  the  following  sections  of  the  paper).  The  coupled  equations  are  given  by 

/ 
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We  found  that  at  any  depth  rjn  the  resulting  iterating  scheme  given  by  Eqs.  (A5)  and  (A6) 


converges  very  fast  with  a  second-order  Taylor  series  expansion  of  the  matrix  exponential  in  Eq. 
(A5).  By  selecting  a  sufficiently  small  At,  we  can  keep  the  number  of  iterations  k  to  2  or  3. 
However,  for  pulses  that  are  short  in  the  temporal  domain,  the  grid  size  must  be  significantly 
increased  to  ensure  that  eigenvalues  of  the  matrix  in  the  exponent  are  less  than  one,  making  the 
Taylor  expansion  possible.  An  alternative  way  is  to  sub-sample  the  grid  on  demand,  i.e.  only  at 
the  high  intensity  areas,  which  are  relatively  small  compared  to  entire  parametric  domain.  At 
every  time  step  rj  l/2  — »  r.+1/2  in  Eq.  (A5),  we  check  the  magnitude  of  the  elements  of  the  matrix 


in  the  exponent.  If  necessary,  we  introduce  additional  M  —  1  time  sub-samples,  such  as 
r;-i/2  =  r°  — »  r1  -»  ...  tm~x  -»  r(.+1/2  =  rM  ,  where  M  is  chosen  such  that  the  elements  of  the 


resulting  refined  matrices  corresponding  to  those  sub-samples  are  small  enough  to  ensure  the 
validity  of  the  Taylor  series  approximation  of  their  matrix  exponentials  [which  are  described 
later  in  Eq.  (A9)].  To  calculate  the  refined  matrices  for  every  such  rm  ?  we  repeat  the  same 
integration  step,  which  was  done  to  derive  Eq.  (A5)  with  the  only  difference  being  that  now  the 
integration  domain  is  given  by  |^r"',  r'"+1  J  .  During  such  integration,  we  also  need  to  estimate  the 

average  intensity  values  at  im=  r"‘  +  I  /  2 A  r'  samples,  where  At'  -  rm+1  -rm  =  At  I M  is  the 
resulting  sub-sampling  time  step.  As  any  such  sample  fm  is  equal  to  r,._1/2  +{m  +  l/2)Ar' ,  we 
estimate  InJ(im)  =  l(r/n,prTm)  by  a  linear  interpolation  of  r,_1/2)  and  /„  .( ri+1/2)  given  by 

Aj  (#■)  =  -1/2 +(*■  +  !/ 2)Ar') 


=  V.  J  (  r,-l/2  )  +  Mj.J  Om/2 )  ■  (A7) 

where  Am  +  jum  =  1 ,  Am  =  (1  -  {m  + 1  /  2)  / M)  .  The  powers  of  the  intensity  values  7“ . ,  a  =  1.3  , 
for  the  half-way  time  grid  samples  r,.  1/2 ,  tm/2  in  Eq.  (A7)  are  interpolated  by 
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Finally,  taking  into  account  Eqs.  (A7)  and  (A8),  for  every  affected  time  sample  the  refined 
version  of  Eq.  (A5)  is  given  by 
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where  M  is  chosen  to  make  the  elements  of  the  matrix  in  the  exponent  of  Eq.  (A5)  smaller  than 
a  given  threshold  i.e.  0  <  s  <  1 ,  such  that 
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The  analysis  of  nonlinear  materials  is  often  guided  by  measurements  of  the  optical 
transmission.  In  most  calculations  the  radial  domain  is  assumed  to  be  constant  but  by  using  both 
the  radial  and  temporal  profiles  of  the  solution,  we  can  now  analyze  and  compare  the 
transmission  plots  of  the  peak  transmittance  Ts ,  defined  by  Eq.  (All),  and  the  integrated 


transmittance  TE,  given  by  Eq.  (A  12)  (with  TF  as  a  cumulative  Fresnel  transmittance  at  the 
interfaces) 
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(A  12) 


The  integrated  value  TE  is  generally  accepted  as  a  key  parameter  quantifying  nonlinear 

materials  because  it  can  be  measured  in  the  laboratory  using  widely  available  thennal  detectors, 
which  average  the  pulse  intensity  over  space  and  time.  However,  as  the  peak  intensity  often 
causes  damage  to  opto-electronic  detectors  and  sensors,  it  is  useful  to  have  a  numerical  check  on 
the  validity  of  the  laboratory  measurements  using  thermal  detectors.  Because  pulse  distortion 
can  occur  in  both  the  temporal  and  radial  domains,  our  calculation  can  also  search  the  radial  and 
temporal  domains  of  the  pulse  to  find  the  maximum  value  of  the  intensity  to  use  in  the 
calculation  of  Ts . 
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Appendix  B.  Description  of  Contribution  to  Each  Energy  Level 


We  define  the  total  intensity  reduction  Qm  for  a  given  grid  index  O)  =  \  n  +  1/2,  /,  /  —  1  /  2  j  as 
the  product  of  individual  reductions  Qs  co  due  to  electronic  levels  s  with  nontrivial  absorption 
cross-sections: 


Q.  =  ||  Q„ 
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sa  ={x|(Ti[v]  +  (T2[x]  +  (T3[x]>0}.  (Bl) 

The  individual  reductions  Qsco  come  from  Eq.  (A6).  They  equal  to  the  corresponding 

exponential  terms  responsible  for  intensity  decrease  at  the  final  iteration,  k  =  K ,  and  are  given 
by: 
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As  we  are  interested  in  relative  contributions,  we  introduce  the  intensity  decay  values  as  follows: 


*„  =  1  -  •  (B3) 

Ps-,o  =  1  -  Q,„  ■  (B4) 

The  total  intensity  decay  values  qco  are  used  to  analyze  what  part  of  the  pulse  mostly  decreases  at 
a  given  depth  value.  Meanwhile,  the  intensity  decay  ps  co  due  to  level  s,  in  turn,  defines  the 
relative  contribution  ps  co  by  the  following  fonnula 


s;o 


(B5) 


Nonlinear  dependency  of  the  relative  contributions  ps  0)  on  the  population  densities  and 


intensities  confirms  that  it  may  be  not  accurate  to  approximate  such  contributions  from  the 
available  values  of  the  population  densities.  In  sections  above  we  first  plot  total  intensity  decays 
q{(0] ,  and  secondly,  we  superimpose  plots  of  absolute  contributions  qs.{co]  derived  from  relative 


contributions  ps  ,0J,  by  scaling  them  to  the  total  intensity  decays: 


?,;{»}  =A,{®}  •?{»}•  (B6) 

By  using  this  feature  in  analyzing  certain  nonlinear  materials,  we  explained  why  our  numerical 
scheme  is  as  good  as  the  other  techniques  in  some  cases,  while  gives  better  solution  in  other 
cases,  and  we  outlined  our  hypothesis  on  why  all  techniques  fail  in  the  case  of  high  input 
energies. 
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Appendix  C  Numerical  Method  in  terms  of  Electric  Field 


84,85 


A  numerical  scheme  to  solve  the  rate  equation  (23)  is  described  here.  An  iteration  fonnula  that 
may  be  used  to  solve  this  equation  numerically  can  be  written  as: 
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where  5Da  =  Ar Da,  a  e  [0..JVJ  .  This  equation  can  be  expressed  in  operator  form  as: 
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with  N(.+1/2  —  [N„+1/2>0>,--1/2 . N)|+1/2j.|j_1/2,...,N)1+1/2^j|._1/2  J  and  ^h+1/2,/,1-1/2  —  ^^Vn+\l2->  Pj^i-Vl)  • 

The  parameter  Yf  [n\e>Da  ~{Q“j  ,■  +  Q“+iji}  may  be  considered  as  a  discretization  of  Y^(//,  r) 

and  can  be  used  in  the  expression  for  the  rate  operator  with  Qn  /(.  =  Q(r/n,pj,ri)  .  A  general 
solution  for  the  propagation  equation  (33)  can  be  expressed  as: 
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A  symmetric  split-step  (Fourier)  technique  can  be  used  to  solve  this  propagation  equation.  A 
value  at  a  subsequent  depth  sample  is  determined  from  the  value  at  a  current  depth  sample  by 
first  applying  a  diffraction  operator  along  A  77  while  neglecting  a  rate  operator,  and  then  applying 
the  rate  operator  to  the  resulting  Q  value  while  neglecting  the  diffraction  operator. 
Mathematically,  this  is  equivalent  to  a  commuting  of  the  rate  and  the  diffraction  operators,  which 
is  expressed  as: 
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Such  operators  may  not  generally  commute,  and  an  error  introduced  by  this  approximation  can 
be  on  the  order  of  A  77 .  To  decrease  the  computational  error,  a  symmetric  version  of  the  split- 
step  procedure  can  be  used,  which  may  be  described  by  the  following  equation: 
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An  expression  for  the  rate-related  operator  is  derived  by  using  the  following  approximation: 
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Substituting  the  approximation  for  the  rate  operator  provided  in  Eq.  (C6)  into  Eq.  (C5)  provides 
a  form  of  the  propagation  equation  which  may  be  expressed  as: 
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A  solution  for  the  diffraction  equation  is  obtained  by  using  a  Crank-Nicholson  scheme.  We  use 
the  following  approximation  to  the  diffraction  operator: 


T#(Ar)*T(p)  =  -V^ 


(C8) 


Applying  a  differential  operator  exp ('P^A 77/ 2) °Q(rj)  is  equivalent  to  solving  the  following 
differential  equation: 


m*’p)  =  ly(p)j6(?,p).  (C9) 

a/7  2 

Substituting  the  expression  for  the  diffraction  operator  provided  in  Eq.  (C8)  yields  the  following 
equation  to  be  solved: 
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where  am  =  1  /  2  x  /  /  4  and  a  time  variable  is  omitted  when  applying  the  diffraction  operator.  A 
forward  difference  approximation  for  the  77  derivative  is  written  as 
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A  variant  of  Crank-Nicholson  scheme  can  be  used  to  approximate  the  Laplacian  operator: 


(Cll) 
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Applying  these  operators  results  in  the  following  iteration  scheme: 
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After  grouping  the  terms  containing  identical  (or  substantially  similar)  Q  samples,  the 
diffraction  equation  can  be  expressed  in  a  discrete  fonn  as: 
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Equations  (Cl 5)  and  (Cl 6)  are  applicable  only  for  “internal”  indices  0<pj<pN, 
0  <  px  <  Pj  <  pN  ,  <  pN  ,  and  boundary  cases  should  be  evaluated  separately. 

A  solution  to  the  differential  equation  that  includes  diffraction  effects  is  obtained  at  a  boundary 
where  p  =  0 .  A  derivative  dQ/dp  can  be  zero  at  p  =  0  in  a  system  having  cylindrical 
symmetry  such  as  a  circular  pulse  footprint.  A  limit  of  the  expression  Q'(p)/  p  is  obtained 
when  p  — »  0  by  using  a  Maclaurin  expansion: 
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Using  Eqs.  (Cl 7)  and  (Cl 8),  the  diffraction  equation  at  p  =  0  is  written  as 
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A  variant  of  Crank-Nicholson  technique  is  applied  to  estimate  the  second  derivative  term  for 

p  — »  0  : 
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A  second-order  approximation  to  the  second  derivative  is  written  as: 
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Using  Eqs.  (C20)  and  (C21),  the  diffraction  equation  at  the  center  of  a  pulse  is  written  in  an 
approximate  form  as: 
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The  derivative  of  (f  can  be  zero  with  respect  to  radial  distance  at  p  =  0  ,  because  the  center  of 
the  pulse  is  an  extreme  value.  Parameter  values  can  be  approximated  at  j= — 1  as  Q" j+1  =  Q[‘"  and 
Q” i  =  Q” .  Therefore  the  diffraction  equation  can  be  rewritten  as  follows: 
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Eq.  (C23)  can  be  used  to  obtain  the  following  general  equation  which  is  applicable  at  a  center  of 
a  pulse  (at  p  -  0  ): 
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A  solution  to  the  diffraction  equation  at  an  outer  boundary  of  a  pulse  at  p  =  pN  ,  can  be  obtained 

using  a  technique  similar  to  that  used  to  find  a  solution  at  the  center  of  the  pulse  described  here. 
A  value  at  a  non-existing  sample  point  p  =  pN  +1  can  be  reconstructed  by  evaluating  a  linearly 

interpolation  at  values  of p  near  the  boundary: 
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This  linear  interpolation  leads  to  a  relationship,  such  as: 
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Eq.  (C27)  can  be  substituted  into  the  general  diffraction  equation,  Eq.  (Cl 5)  to  provide  a 
diffraction  equation  having  a  form 
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Eq.  (C28)  can  be  expressed  in  a  simpler  form  as: 
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Eqs.  (C15),  (C24)  and  (C29)  can  be  combined  to  provide  a  propagation  procedure  that  includes 
diffraction  effects,  which  can  be  calculated  as  a  solution  of  the  following  linear  system  of 
equations  with  a  tri-diagonal  matrix: 


where 


and 


UQ"+1  =  VQ" , 


q* =[&'- 
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(C31) 
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The  matrix  V  can  be  expressed  in  a  form  analogous  to  that  of  the  matrix  U  in  Eq.  (C33).  A 
conventional  technique  applicable  to  tri-diagonal  matrix  systems  can  be  used  to  solve  Eq.  (C31) 
such  as,  e.g.,  an  LEI  decomposition  using  forward-  and  back-  substitution  for  a  tri-diagonal 
system.  A  solution  to  a  linear  equation  such  as  that  in  Eq.  (C3 1 )  is  stable,  and  a  zero  pivoting  is 
not  required  for  a  matrix  which  has  a  diagonal  dominance  at  its  rows.  A  diagonal  dominance 
condition  can  be  expressed  for  a  tri-diagonal  matrix  by  the  relationship  as  follows: 


KI>Kil+|<i|.  (C34) 

Elsing  the  approximation  pj  =  jAp ,  the  condition  provided  in  Eq.  (C34)  generally  holds  for  the 
central  and  internal  sample  points: 
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(C35) 


However,  Eq.  (C35)  uses  a  constraint  on  incremental  delta  values  used  because  of  a  condition  at 
a  boundary  sample  point: 
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which  may  be  true  if: 
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which  can  provide  a  constraint: 
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A  solution  to  the  propagation  equation  provided  in  Eq.  (C7)  uses  a  specification  of  an  operator 
Qp(rj  +  A 77,  r)  ,  which  depends  on  an  unknown  //  +  A //  sample  point.  An  iteration  procedure 

should  be  applied,  where  this  sample  point  at  a  current  iteration  is  drawn  from  results  of  previous 
iteraiotns.  A  detennination  of  a  subsequent  sample  point  value  for  Q  can  be  written  as: 


Q"+1  <- 


TV"+l/2«  y  1 V«+1/2«  ] 

1/1 1+1/2  Y  M/2 
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According  to  Eqs.  (C1)-(C3),  an  iteration  step  for  the  rate  equation  has  a  form: 
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^wM  =  T0+XTfM, 


(C41) 


S,}- 
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Each  iteration  involves  applying  three  operations:  applying  the  diffraction  operator  at  k  =  0 
using  Eq.  (C43),  applying  polarization  term  using  Eq.  (C44)  and  applying  the  same  diffraction 
effect  using  Eq.  (C45).  Eqs.  (C43)-(C45)  are  written  as 


Qmid  =  U'‘VQ"  , 


Q;:+1(A+1)=exp(-A70jnexp(A70f  [n  +  l/2,/-l/2]0f  [n])Q: 


(C43) 

(C44) 
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where  the  following  discrete  forms  of  propagation  operators  0^(77,  r)  and  0^(7/,  r)  are  used: 
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of'  [n  +  \l2,i-\l2]  =  <sp- 
^o{k]  r  i  Q 

H=- 

and 

A  computational  procedure  for  determining  the  propagation  of  an  electromagnetic  pulse  in  a 
nonlinear  absorbing  material  which  accounts  for  diffraction  effects  are  described  using  the 
following  equations.  Based  on  Eqs.  (C1)-(C2),  or  (C40)-(C42),  a  rate  equation  are  used  which 
have  the  following  form: 
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i0  =  A zT0D0 ,  and  ia  =  A  r  ro  ( 2nc0g02 )  .  (C50) 

anco0  v  7 

Using  the  results  provided  in  Eqs.  (C43)-(C48)  a  mathematical  expression  to  describe 
propagation  of  a  pulse  can  be  written  in  the  following  form: 
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By  introducing  a  parameter  %  = 


a,„A7 

A P1 


the  nonzero  elements  of  internal  rows  of  matrices  U  and 


V  ,  i.e.  u'j  and  v'.  such  that  i  e  Nr }  and  j  e  {/ - 1, . . . , i  + 1} ,  can  be  written  as: 
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vU=~uU  ’ 

(C58) 

vJj=\-Z  =  2-uJj  ; 

(C59) 
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By  applying  Eq.  (C25),  the  nonzero  elements  of  the  first  rows  of  matrices  U  and  V  ,  i.e.  and 
v'j  such  that  i  -  0  and  j  e  {0,1} ,  can  be  written  as: 


n0°=  l  +  2£  ;  (C61) 

Mj°  =-2£  ;  (C62) 

v00=l-2£  =  2-n0°  ;  (C63) 

and 

v,°  =  2£  =  .  (C64) 

By  applying  Eq.  (C30),  the  nonzero  elements  of  the  last  rows  of  matrices  U  and  V  ,  i.e.  u'  and 
v'j  such  that  i  =  Nr  and  j  e  {Nr  - 1,  AT.}  ,  can  be  written  as: 
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